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Summary 


Variational calculus methods are applied to the problem of torsional 
oscillations of elastic spheres of variable densities and elasticities and 
the corresponding plane Love wave problem for vertically heterogeneous 
media. ‘The methods work well and give the phase velocity curve shown 
in Figure 1 tor the model Earth, the internal constitution of which was 
interred by K. E. Bullen. “Vhe corresponding displacement distribu- 
tions are shown in Figure 2 


1. Introduction 


‘Torsional oscillations of elastic spheres of variable densities and elasticities 
have been studied by several authors. Recently, Matsumoto & Sato (1934) have 
studied the torsional oscillations of a sphere consisting ot the mantle and core, each 
ot which is homogeneous. N. Jobert (1956) has studied the torsional oscillation 
of spherical harmonic degree 2 type of the Earth, the internal constitution of which 
was inferred by K. E. Bullen. Making use of the calculus of variations, she has 
got the period of 43°54min for the Earth thus assumed. All of these authors are 
concerned with the oscillations of rather lower modes. As was pointed out by 
J. TL. Jeans as early as in 1923, by making the wave number larger and larger in 
this kind of problem, we arrive at the corresponding surface wave problem. ‘hus 
the study of torsional oscillation of the Earth mav throw some light on the nature 
of the G wave. 


2. Variational method 


From the mathematical point of view, the problem of torsional oscillations of 
elastic spheres of variable densities and elasticities is reduced to solving the 
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Summary 


Variational calculus methods are applied to the problem of torsional 
oscillations of elastic spheres of variable densities and elasticities and 
the corresponding plane Love wave problem for vertically heterogeneous 
media. The methods work well and give the phase velocity curve shown 
in Figure 1 for the model Earth, the internal constitution of which was 
inferred by K. E. Bullen. The corresponding displacement distribu- 
tions are shown in Figure 2. 


1. Introduction 


Torsional oscillations of elastic spheres of variable densities and elasticities 
have been studied by several authors. Recently, Matsumoto & Sato (1954) have 
studied the torsional oscillations of a sphere consisting of the mantle and core, each 
of which is homogeneous. N. Jobert (1956) has studied the torsional oscillation 
of spherical harmonic degree 2 type of the Earth, the internal constitution of which 
was inferred by K. E. Bullen. Making use of the calculus of variations, she has 
got the period of 43°54 min for the Earth thus assumed. All of these authors are 
concerned with the oscillations of rather lower modes. As was pointed out by 
J. H. Jeans as early as in 1923, by making the wave number larger and larger in 
this kind of problem, we arrive at the corresponding surface wave problem. Thus 
the study of torsional oscillation of the Earth may throw some light on the nature 
of the G wave. 


2. Variational method 


From the mathematical point of view, the problem of torsional oscillations of 
elastic spheres of variable densities and elasticities is reduced to solving the 
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differential equation 


under the boundary 


(2.2) 


i!) and / are continuous at = &p, say. (2.3) 


In (2.1), p and p are density and rigidity, each being assumed to be a function of 
&, n the degree of spherical harmonic, p the frequency, e*?¢ being understood 
throughout the equations. ¢ is the distance r from the centre of the sphere divided 
by the radius of the sphere a. £9 in (2.3) is € corresponding to some surface 
of discontinuity within the sphere, if one exists. The first and second conditions 
in (2.3) correspond to the continuities of displacement and stress, respectively, 


whereas (2.~ shows the vanishing of stress at the free surface r = a. (2.1)-(2.3) 
will be obti..ued by putting 


r=at, (2.4) 


in (102) and (104) in a previous paper (Takeuchi 1950). 
We can easily find the variational equation equivalent to (2.1)-(2.3), which, 
after a little transformation, is reduced to 


1 
| de 
0 


As the trial function in (2.5), we take 


where the A» are some constants to be determined and each [(m, £) or simply 
K(m) is chosen so as to satisfy the conditions (2.2) and (2.3). We understand also 
that the variation process is applied only to the undetermined constants Am. 
In this case, (2.5) are reduced to m equations of the following type 


for the assumed /(m) of number m in all. In fact, putting (2.6) into (2.7), we get m 


and 
| 
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simultaneous equations among m unknown Am. Putting the determinant formed 
by the coefficients equal to zero will give the value of p required. 

In our calculations, we take into account the existence of the fluid core and 
assume no surface of discontinuity within the mantle. Thus the condition (2.3) 
becomes 

dl ani =o at é 6370 — 2900 
dé 6370 


As the trial function /(m) which satisfies (2.2) and (2.8), we take 


m—1 347 mi 347 

where m is some positive integer. Since the boundary conditions are independent 
of the wave number 1, /(m) in (2.9) may be a trial function for any n. This will make 
our calculation quite short. In our practical calculations, we use only three func- 
tions [(2), (6) and (20). The reason why we use discrete values of m as above is 
as follows. As is seen after a little numerical calculation, /(3) and /(4), for example, 
will give almost the same distribution (mode) as (2). Thus, if we use /(2), (3) and 
(4) in our calculation, for example, we lose almost all significant figures and get 
no reliable final result. In the calculation of [(2), (6) and U(20) also, we must take 
care not to lose significant figures. This we can do by keeping fractional expres- 
sions for numbers as late as possible. 


= 0°544741. (2.8) 


(2.9) 


3- Trial results 

Using the method described in the last section, we shall at first check the 
results obtained by Matsumoto and Sato. The calculation proceeds as follows. 
We at first calculate the values of Uz), (6), (20) and their derivatives at € cor- 
responding to the depths of 0, 200, 400, 700, 1 000, 1 400, 1 800, 2350 and 2g00km 
in the Earth. Next, using Simpson’s method of numerical integration, we calculate 
the integrals in (2.7). In the present case, p and yu are assumed to be constant, 
so they go out of the integral. The results of calculations are expressed as follows. 


22 2 Ag+20 2 Ax = 0 
26 Ao+6 6 Ag+20 6 Ax = 0 (3-1 


220 Ao+620 Ag +20 20 Ag = 0 


where the coefficient 2 2, for example, is 


2 2 = 47°8392 P—159°9777N + 319°6554 
N=n(n+1) P= » Veo = J(®) (3-2) 
po 


po and po being p and p at the surface. 2 2, together with other coefficients, is 
shown in a self-explanatory way in Table 1. The coefficients are symmetrical with 


respect to the diagonal. So we have, for example, 26=62. Putting the deter- 
minant formed by these coefficients equal to zero, we have the eigen value P for 


Coy 


given N. We shall call the results obtained by using 2 2; 2 2, 2 6, 6 6; 2 2, 


92 H. Takeuchi 


2 6, 6 6, 2 20, 6 20, 2020 as the first, second and third approximations, respec- 
tively. The first, second and third approximations of P for n = 2, 3, 4 and oo 
are shown in Table 2. In Table 2, are shown also the exact P obtained by Matsu- 
moto and Sato. We see that our second approximations give values of P almost 


Table 1 
P N 


47°8392 —159°9777 319°9554 

—1°013227 1°703042 —3°406082 

—0'5048223 08250114 —1°650023 
00377500 —0'0693672 —2:°188851 
0019604738 —0°02950595 —1°317838 
0:011661966 —0°01950237 —1°621529 


identical to theirs. Our second approximations give one more P. They are 
P = 8-758, 9°372 and 10°32 for nm = 2, 3 and 4, which agree quite well with the 
exact values P = 8-18, 8-9 and g-g obtained by Matsumoto and SAato. By the 
reason which will be given later, we must have P = 1-0 n as n tends to infinity. 
Our second and third approximations give P = 1-075 n and P = 1-032 n as n 
tends to infinity, which are not so good. The reason why we get rather poor result 
in this case is simply that ((m) with m = 2, 6 and 20 are not good approximations 
to the mode for this case. We are sure to get better approximations by using /(m) 
with larger m, but we shall not try to do so here. As is well known, approximate 
values by using the variational calculus here are always larger than the true values. 
We see in the above this is always the case. 


Table 2 


pa 27a 
P= —= 


Vio 
Ist 2nd 3rd Exact 


2°423 2°413 2°415 
3°831 3°795 3°79 
4°99 
1-212" 1032" = (10 n) 


4- Results for Bullen’s Earth model 


Having seen the efficiency of our method in this kind of calculation, we shall 
now calculate P for the model Earth, the internal constitution of which is inferred 
by K. E. Bullen. The numerical data are taken from Bullard’s text book (1954) 
and are shown in Table 3. 

In Table 3, surface values are obtained by extrapolation from those in the deeper 
part. Thus we do not take into account the existence of the Earth’s crust. This 
may be allowed in this kind of calculation. The calculation proceeds as in the last 


l 
a2 
26 
2 20 
6 6 
6 20 
20 20 
2 
3 
4 
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Table 3 


Properties of Earth model 


Depth = a—r p 
km g/cm 10!2 dyne/cm 


° 
200 0-74 
400 0-89 
700 1°52 

1 000 ‘ 1°89 

1 400 . 2°15 

1 800 2°39 

2350 2°73 

2.900 3°03 


section except that we have now variable density and elasticity. Corresponding to 
Tables 1 and 2, we have now Tables 4 and 5. 


Table 4 
P N l 


36-3870 —53°4030 106-8061 
—0°874861 0982992 —1°965984 
—0°452904 0°5080651 —1:°016130 

00330281 —0:0382040 —0°755642 
©°01792547 —0°01974089 —0°6230995 
0°010914097 —0'01177460 —0°°974379 


Table 5 


Convergence of approximations to P 


3rd 


7°285 
12°92 
22°29 
39°50 

1127" 


In the present case, we have no exact value to compare with. However, as the 
successive approximate values converge quite well, we are sure to get reliable 
approximations for P, especially for smaller n. In the present case also, we must 
get P = 1-0 m as n tends to infinity. Our third approximation gives P = 1-127n 
as n tends to infinity, which is not so good. The reason why we get poor results 


22 
ory 
26 
2 20 
> 
66 
6 20 
om 
20 20 
pa 27a 
P 
TV. 
n Ist 2nd || 
2 3°658 3°601 
4 7°758 7°287 
8 15°30 12°97 
16 30°05 23°34 
32 59°37 44°52 
a 1‘829n 1°353” 
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here again is the same as at the end of the last section. This will be understood as 


follows. Putting the coefficient 20 20 for n = 32, for example, equal to zero, we 
get P = 39°74, which is very close to the third approximation P = 39-50 shown in 
Table 5. This result shows that ((m) with m = 20 is a good approximation to the 
mode for n = 32 and we shall have a better result by using this one in the first 
approximation and some other /(m) with m larger than 20 in the second and third 
approximations. However, we shall not try to do so here. We know in the next 
section that even the results for n = 16 and 32 in Table 5 are good enough for 
our present purpose. Putting a = 6370km and V5.9 = 4:26km/s, which is obtained 
from yu and p at the surface, into the third approximations in Table 5, we get the 
periods T shown in Table 6. We see that the period for n = 2 is 43-4min, which 
agrees well with 43:54min obtained by N. Jobert. The asymptotic expression 
for the Legendre function P(@) shows that the wave length A for large n is equal 
to 27a divided by n+ }. Dividing this wave length by the period in Table 6, we get 


Table 6 
Wavelength, 4, Phase, C, and Period, T 

n A 

km km/s min 
2 16 000 6°140 43°43 
4 8888 6-903 21°48 
8 4706 6-481 
16 2424 5°761 7°02 
32 1231 5°181 3°96 


the corresponding phase velocity C. Although the above procedure is valid only 
for large n, we show in Table 6 the values of A and C obtained in this way. The 
reason why we must have P = 1-on as n tends to infinity is now clear. Since 
P/(n+4) is equal to C/Vs.9 the above requirement is another way of stating the 
well known fact that the phase velocity of Love waves for very short wave length is 
equal to the shear wave velocity at the surface of the Earth. 

The results in Table 6, together with some of the results which will be obtained 
in the next section, are shown in Figure 1. From this figure, we see that the phase 
velocity C has a maximum of about 6-9km/s for the period T of about 1 300s 
or for the wave length A of about goookm. That C becomes smaller for wave 
length longer than this critical one shows the existence of some lower (shear wave) 
velocity layer in the deeper depth. In the present case, this lower velocity layer is 
nothing but the Earth core. The core boundary is at 2gookm depth. The scale 
ratio, wave length/depth, is about 3, a plausible one. In the following way, we get 
the scale ratio of this order of magnitude for other wave lengths. The phase velocity 
C for n = 8 (or A = 4706km) is 6-481 km/s. The depth at which the shear wave 
velocity takes this value is at about 1 200km within the Earth. The scale ratio in 
this case is about 4. Similarly, for m = 16 and 32, we have the scale ratios of about 
4 and 3. This kind of scale ratio consideration is sometimes very useful, especially 
when we want to get a rough idea of the dispersion curve for complicated 
structures. 

In this connection, it will be interesting to see the displacement distribution 
with the depth, or vibration mode for each n. This we can do by solving (3.1) 
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and putting the results thus obtained into (2.6) and finally into 


4) 
uy = O, é 20 (4-1) 


u,, ug and ug in (4.1) show (r, 8, 6) components of displacement in spherical co- 
ordinates. Since we are treating here the free oscillation problem, we cannot 


Phase Velocity: km/s 


500 1000 1500 2500 
Period : sec. 


Fic. 1.—Theoretical phase velocity curve of the mantle Love wave. 
expect to get an absolute value for the displacement. Thus we show in Figure 2 


the distribution of //£ against depth, taking its value at the surface as unity. We 
see from this figure that the shorter the wave length, the shallower the depth to 


10 


Displacement 


Depth: km. 


Fic. 2.—Displacement distributions within the Earth. n is the degree 
of the corresponding spherical harmonic. 


which the corresponding displacement penetrates. This, of course, must be the 
case. For very short wave lengths, we have surface waves in the usual meaning. 


| 
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For n = 8, 16 and 32 in Figure 2, the displacements almost vanish within the 
mantle. Thus for these waves or oscillations, the existence of the core of the Earth 
has almost no influence. For n = 4, the displacement at the core boundary is about 
} of that at the surface, which is not negligible. ‘Thus in this type of oscillation we 
can expect some effect of the existence of the Earth’s core. This is exactly what we 
had in the above. 

To find values of C not calculated above, Figure 1 may be used. We can find 
C for n = 6, for example, as follows. First, we have the wave length A = 6154km. 
We try next to find C and T which are on the smooth curve in Figure 1 and the 
product of which is equal to 6154km. Thus we get T = goo, C = 6-75 approxi- 
mately. 


5. Variational method for Love waves 


In order to make connection with the usual surface wave problem, we shall 
study in this section the plane Love wave problem in vertically heterogeneous 
media. There are many studies on this problem. References may be made to the 
text book by Ewing & others (1957). Putting 


vo etk(Ct—2z) (5.1) 


in their equation (7.69), we get the following fundamental equations to be studied. 


dy dp do 
— + uly = 0. : 
+ + = 0 (5.2) 


The boundary conditions are the vanishing of stress at the free surface z = 0 
and continuities of displacement and stress at some surface of discontinuity if it 
exists. These conditions are satisfied by putting 


dv 
—=0 at z=0 (5-3) 
dz 
and 
dv 
v and continuous (5.4) 


at the surface of discontinuity. We can easily find the variational equation (Jeffreys 
1934) corresponding to (5.1)~(5.4), which, after a little transformation, is reduced to 


+00 +00 
1 dv div 
C2vbvdz— ov + — — =o. 


We consider the half-space in the + 2 region and assume that there is no surface 
of discontinuity there. We proceed as in Sections 2 and 3. We first put 


= Am o(m, ke) (5.6) 


where Am are some undetermined constants, and as v(m, 2), each satisfying (5.3), 
we shall take 


kz) = e-mk2 — o-5 e-2mkez, (5-7) 


Variations in A» will give equations similar to (3.1). The coefficient of A,’ 
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in the mth equation may be put in the form 


du(m', kz) do(m, kz) 
kz) — \0(m’, kz)o(m, kz) + (5.8) 


The reason why we take kz as the independent variable rather than z is that 
we can then carry out almost all calculations independently of the wave number k. 
The consideration of k comes only at the last step when we form the products 
including p and pw. The distributions of p and yw in the kz space are different 
for different k. In our practical numerical calculus, we use m = 0°5, 1-0 and 
sometimes m = 0°25, 2:0. 

In order to check the efficiency of our method, we first apply it to the case 
studied by Wilson (1940, see also the text book by Ewing & others, p. 214). That 
is the typical Love wave problem, a superficial layer on a semi-infinite medium. 
The rigidity in the lower medium is 1-8 times that in the upper layer. The shear 
wave velocities in the upper and lower layers are 3-6 and 4-6km/s, respectively. 
As there is a surface of discontinuity in the region considered, our method described 
above seems not to be applicable to the present case. But by replacing the discon- 
tinuity surface by a thin continuous layer in which rigidity and shear wave velocity 
change rapidly, we may apply our method here. Taking v(0°5, kz) as a trial 
function, and carrying out the integration, we get the results in Table 7. 


Table 7 


| (—) d(kz) 
kn 


0°00 


1°72 0°023215956 0°0051717190 


2°84 
4°27 
6°86 
10°02 
16°24 
fore) 


0°086799923 
0°16283787 
0°25141526 
©°30955991 
0°36368226 
0°45833333 


0'016788801 
0°027316495 
0°035677655 
0039073316 
0°040879418 
0°041666667 


Table 7 shows that our approximate method gives almost identical results with 
Wilson’s except for kH - o and kH -> oo, thus justifying the replacement of the 
surface of discontinuity by a continuous layer. In order to remedy rather poor 
results for kH - 0 and -> 00, we have only to replace our trial function v(0°5, k=) 
by v(0°25, kz). In fact by using v(0-25, kz), we get for RH — o and — o0, C = 3°64 
and 4:65 km/s, respectively, which agree well with the exact results C = 3-60 and 
4°60 km/s. Thus we may expect very good agreement in the second approximation 
by using v(0-5, kz) and v(0-25, kz). We shall, however, not try to do so here. 
For our present purpose, it is sufficient to know that our method works well. 


6. Results for mantle waves 


Having thus checked the efficiency of our method, we shall now apply the 
method to the plane boundary problem corresponding to the case discussed in 


on 
| v2d(kz) 
kH als Our C Exact C 
km/s km/s 
3°76 3°60 
3°83 3°80 
4°02 4°00 
4°22 4°20 
4°42 4°40 
4°54 4°5° 
4°64 4°56 
4°81 4°60 
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Section 4. The density and rigidity distributions are shown in Table 3. We restrict 
our study to the case when there is no need to know the distributions within the 
core. The computation proceeds as follows. We first calculate the values of 
v(0°5, kz), v(1-0, kz), v(2-0, kz) and their derivatives (with respect to kz) for several 
values of kz, and then form the products like v(o-5, kz) v(1-0, kz) and 


dv k dv he) 
*O, RZ). 
dha? 


Substituting values of k will give the corresponding distributions of p and yp for 
each kz. Forming products like pv(o-5, kz) v(1°0, kz), we carry out the numerical 
integration. The coefficients, which have meanings similar to those in Table 4, 
for 1/k = 200, 100 and 50km are shown in Table 8. The wave lengths for these 
waves are 1257, 628 and 314km respectively. 


Table 8 
(C/ I (C/ V,.0)? I (C/ I 
O°5 O'5 1°495349 —2°164886 1425821 —1°764019 1400108 = —1°625459 


O°5 0°9490493 —1°327204 09246861 —1°175540 0°9143637  —1°114840 
1°O 0°7157905 —1°126688 0°7031800 —1°'031734 0°6972855 —0:9883948 


1/k 200 km 100 5° 


In the similar way as in Sections 3 and 4, we have the first and second approxima- 
tions of C/V, 9 for each k as shown in Table g. Since they show fairly good conver- 
gence, multiplying the second approximation by V, 9 = 4:26km/s, we may get 
the value of C for each k as shown also in Table 9. These results are shown 


Table 9 
1/k C/V;.9 C A T 
Ist 2nd 
200 km 1'203 I°199 km/s 1257km 246s 
100 4°73 628 133 
50 1069 4°55 314 69 


in Figure 1 by black circles. The smooth curve in this part of Figure 1 is drawn 
by connecting the point corresponding to n = 16 and the point (T = 0, C = Vyo 
= 4:26km/s). The merging of the results obtained in this Section and in Section 
4 into each other will assure us of the validity of our methods of study. 

Figure 1 shows that the phase velocity C is almost linear against the period T 
for T less than 400s. This result agrees well with the study of G waves by Sato 
(1958). In this connection, it will be interesting to work out the corresponding 
group velocity U. Assuming that C is linear against T between T = 0 and 
4008, we get 

dC = 4°26 


aT = 0°0036km/s?. (6.1) 
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The formula 
U= (6.2) 
1+(T/C\(dC/dT) 

then gives U = 4:29, 4°35, 4°44 and 4°55 km/s for T = 100, 200, 300 and 400s. 
These are about the values of group velocities usually observed in this period inter- 
val. We must take care, however, not to attach much importance to these figures. 
The calculation has not taken into account the existence of the Earth’s crust. 
Furthermore, V,, = 4:26km/s in the Jeffreys and Bullen model is a little far 
from the truth. Recent studies in explosion seismology indicate that V, 9 is about 
4°7km/s. In the present paper, we are concerned with rather large-scale charac- 
teristics of the Earth, so it is sufficient for our present purpose to get figures roughly 
in accord with the observational data. The approximate method in this and the 
last Section, however, is one of the most effective in discussing minor structures in 
the upper mantle and the crust. So after getting enough data, we hope to study 
this interesting problem in the near future. Taking care not to apply rashly the 
usual wave concepts to the vibrations of the longer periods, we calculate the values 
of U for T = 600, 800, 1000 and 1200s. Use is made of Figure 1 and (6.2). 
The results are U = 5-15, 5°71, 6:16 and 6-54 km/s for the above periods, respec- 
tively. No flattening of the dispersion curve is seen here as is found in the mantle 
Rayleigh waves (Ewing & Press 1954, Benioff & Press 1958). The flattening of 
the dispersion curve of mantle Rayleigh waves around the period of 500s is 
considered to be due to the existence of the core. By the present study, this in- 
ference is a little doubtful. The wave of period 500s is of too small scale to be 
influenced by the existence of the Earth’s core. As is said in Section 4, it is only 
atn = 4 or T = 1300s that the wave or vibration is influenced by the existence of 
the core. It may be said against this argument that the situation is different in the 
case of the mantle Rayleigh waves from the mantle Love waves studied here. In 
order to make this point clear, it will be best to study the statical deformation of 
the Earth as a whole by normal traction at the surface. This is a deformation of 
the dilatational type, and as the statical deformation penetrates deeper than the 
corresponding dynamical deformation, the study will be one of the crucial tests 
for the above problem. This kind of study will be made in a forthcoming paper. 


7. Conclusions and acknowledgments 


Periods of oscillations and corresponding phase velocities have been calculated 
for the torsional oscillations of the Earth, the internal constitution of which was 
inferred by K. E. Bullen. The fundamental period calculated is 43-4 min agreeing 
closely with that obtained by N. Jobert. By making the wave number larger, we 
get the phase velocity curve for the G wave. The calculated results agree roughly 
with the observations. The variational calculus method used in the present study 
is very effective and makes the problems considered here matters of computational 
interest. Modern digital computers will be of help in making the computations. 
In connection with the present study, several problems suggest themselves. 
Among them are the problem of statical deformations of the Earth as a whole 
by normal traction at the surface, the study of structures in the crust and in the 
upper mantle and the application of the variational calculus method to the study 
of dilatational oscillations of the Earth and the corresponding mantle Rayleigh 
waves. 
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Summary 

A numerical method for calculating the seismogram of any seismo- 
graph for an arbitrary continuous ground disturbance is used to investi- 
gate the effect of instruments with various physical constants on both 
theoretical and actual impulsive signals. This calculation solves the 
equations of motion of the pendulum-galvanometer system by a 
fourth-order Runge—Kutta process. Variation of instrumental constants 
produces conspicuously different seismograms, in several cases of 
practical interest. 

In order to resolve the apparent discrepancy between the observa- 
tions of Saté, and those of Benioff and Press for the group velocity of 
long period Love waves, certain computed seismograms of the long 
period Love wave arrival, G1, from the New Guinea shock of 1938 
February 1, were compared with the linear strain recording of this 
signal. The results of this study confirm the observations of Saté, 
that the group velocity of Love waves is essentially flat in the long 
period range to several hundred seconds. The distortion introduced by 
instruments whose magnification changes rapidly with period is great 


enough to mask any small dispersion that might be associated with these 
pulse-like signals. 


1. Introduction 


In order to determine precisely the effects of seismographs with various 
physical constants on the pulse-like signals from earthquakes, a systematic investi- 
gation has been made of the modification of certain theoretical pulses, as well as 
signals from earthquakes, using an IBM 650 electronic data processing machine. 
Modifications of seismic ground motion, which are produced by the recording 
instrument, have been discussed previously in many studies of near earthquakes, 
and of teleseisms, for both body and surface waves.* However, the method 
presented here is more powerful, and does not require the use of the arbitrary 
integration constants which are needed to reduce a seismic trace to ground motion. 

Although instrumental distortion affects the entire seismic signal, a quantitative 
measure of this effect is especially important at present in the case of surface 

* See Akima & Nagamune (1956), Birkenhauer (1945), Bryan (1936), Byerley (1942), Coulomb 
(1952), Dyk (1935), Ewing & Press (1952), McComb & Ruge (1937), McComb & Wenner (1936), 


Neumann & others (1943), Ruge (1936), Sat6é (1958), Sharpe (1935), Tsuboi (1931), Walter 
(1945), Wilson (1932), Wilson (1940). 
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waves, because for these waves signal character as well as arrival time is used to 
study the physical properties of the Earth. In one study of well-dispersed Rayleigh 
wave trains, it was shown (Ewing & Press 1952), after comparing data from 
seismograms with different physical constants, that it is neither necessary to 
integrate the seismograms to true ground motion nor to eliminate instrumental 
phase shifts for the periods studied and the precision required. 

There are, however, surface wave signals with small dispersion. A particularly 
striking example of a pulse-like teleseismic signal is the long-period Love wave, 
named the G wave by Byerly (1926) after Gutenberg, for his pioneer work on this 
phase (Sieberg & Gutenberg 1924). His studies, as well as most other investiga- 
tions of this wave* indicated that it shows little dispersion, since its pulse-like 
character is preserved at great distances. These studies demonstrated that the 
Love wave group velocity curve is essentially flat in the long period range to several 
hundred seconds, with a group velocity near 4-48km/s. For a summary of disper- 
sion data on Love waves of shorter periods, see Ewing, Jardetzky & Press (1957). 

In a recent study, Benioff and Press (1958) suggested, on the basis of seismo- 
grams from a new instrument, that the previous studies of the G wave had resulted 
in poor group velocity determinations, and that the Love wave group velocity 
curve has a true maximum, with a decrease in velocity for periods greater than 
about two minutes. The more pronounced dispersion of mantle Rayleigh waves 
which show a maximum for periods near one minute and a decrease in velocity to a 
minimum for periods near four minutes, has permitted determination of the group 
velocity and associated maximum and minimum to be made from the seismograms 
of many instruments (Ewing & others 1957, Ewing & Press 1954(a), (b)). In 
addition, the group velocities obtained from studies of the direct wave agreed 
with determinations based on multiple passages of these waves around the earth, 
although the wave train had been strongly altered by its conspicuous dispersion. 
The agreement of these determinations shows that effect of the instruments 
on mantle Rayleigh waves is small compared to that of the dispersion. 

The essentially constant G wave group velocities were determined by Satd 
(1958), exclusively by the measurement of differences in phase which were obtained 
by Fourier analysis of the arrivals G3, G1, and G4, G2. In this manner, the effects 
of the instrument and of the spectrum of the source are completely removed. 

The group velocities of long-period surface-wave signals provide a means of 
studying the mantle which is independent of body wave studies, and it is therefore 
quite important to get reliable surface wave group velocity data, as the derived 
mantle structures are sensitive to changes in the observed dispersion curves. 

After an investigation of the effect of seismographs on theoretical signals, 
calculations were made of the seismograms which would be produced by various 
instruments upon reception of the G1 signal from the New Guinea earthquake of 
1938 February 1, which had been used by Saté6 for his group velocity determination. 

Our calculations have shown that the instruments significantly influence the 
shape of pulse-like signals (see Figures 1 and 3) in cases of practical interest. Such 
changes in wave shape will cause measurable changes in the observed group 
velocity for signals with slight dispersion unless instrumental effects are removed. 
The modification of the G1 signal by the four pendulum-galvanometer systems 
described in Figure 4, which is similar to the distortion of the test signals in Figure 
3, indicates that the effect of the recording system on these pulse-like signals is 


* See Coulomb (1952), Di Filippo (1950), Ewing & Press (1956), Imamura (1937), Sat6 (1958), 
Sieberg & Gutenberg (1924), Wilson (1940). 
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quite important. Since this is the same G wave signal which, after elimination of 
the effects of the recording system and of the spectrum of the source, did not 
produce a maximum in the group velocity (Sat6é 1958), the earlier determinations 
of an essentially flat Love wave group velocity curve for long periods up to several 
hundred seconds are confirmed. This Love wave group velocity curve is useful 
for testing structures proposed for the upper mantle. An extensive programme of 
study of the upper mantle, using surface waves and body waves is nearing com- 
pletion at the Lamont Geological Observatory, and the results will soon be pub- 
lished (Landisman, Sat6é & Ewing). 

The present investigation of seismograph distortion of pulse-like ground 
motions indicates that, for best results, dispersion curves for pulse-like signals 
should be constructed from the seismograms of instruments whose magnification 
changes very gradually with period. In addition, this study emphasizes the fact 
that the transient response of instruments whose magnification changes rapidly 
with period can be important in studies of wave propagation. 


2. Method 


The transient response of various types of seismographs to pulse-like signals 
will be considered. The equations of motion of these systems, with suitable initial 
conditions, provide a complete description of their behaviour. 

One type of system is the linear strain as recorded by a galvanometer by means 
of a velocity pickup. The response, R, of a directly recorded strain rod, for har- 
monic components much longer than the rod, has been shown by Benioff (1935, 
1955) to be . 


= — sina cosa (1) 


for transverse waves, and 
oz 
-— ats — 
L cos* (2) 
for longitudinal waves, 


where R = displacement of the free end of the horizontal strain rod relative to 
the earth beneath the free end, 
z = horizontal particle displacement, 
L = length of strain rod, 
c = horizontal phase velocity, 


a = angle between strain rod and the azimuthal direction of approach 
of the waves. 


Figure 1 in the study of Benioff & Press (1958) shows that in practical cases, 
changes of phase velocity of the order of 1km/s produce negligible changes of 
magnification. When phase velocity is treated as a constant and the time derivative 
of R, as obtained from a velocity pickup, is recorded on a galvanometer, the magni- 
fication is similar to that of a simple pendulum with the same period and damping 
as the galvanometer. The pendulum obeys the well-known equation 


¥+ 2e9% + wo2x = —2, (3) 


where x = pendulum displacement with respect to the ground, 
z = absolute ground displacement, 
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waves, because for these waves signal character as well as arrival time is used to 
study the physical properties of the Earth. In one study of well-dispersed Rayleigh 
wave trains, it was shown (Ewing & Press 1952), after comparing data from 
seismograms with different physical constants, that it is neither necessary to 
integrate the seismograms to true ground motion nor to eliminate instrumental 
phase shifts for the periods studied and the precision required. 

There are, however, surface wave signals with small dispersion. A particularly 
striking example of a pulse-like teleseismic signal is the long-period Love wave, 
named the G wave by Byerly (1926) after Gutenberg, for his pioneer work on this 
phase (Sieberg & Gutenberg 1924). His studies, as well as most other investiga- 
tions of this wave* indicated that it shows little dispersion, since its pulse-like 
character is preserved at great distances. ‘These studies demonstrated that the 
Love wave group velocity curve is essentially flat in the long period range to several 
hundred seconds, with a group velocity near 4-48km/s. For a summary of disper- 
sion data on Love waves of shorter periods, see Ewing, Jardetzky & Press (1957). 

In a recent study, Benioff and Press (1958) suggested, on the basis of seismo- 
grams from a new instrument, that the previous studies of the G wave had resulted 
in poor group velocity determinations, and that the Love wave group velocity 
curve has a true maximum, with a decrease in velocity for periods greater than 
about two minutes. The more pronounced dispersion of mantle Rayleigh waves 
which show a maximum for periods near one minute and a decrease in velocity to a 
minimum for periods near four minutes, has permitted determination of the group 
velocity and associated maximum and minimum to be made from the seismograms 
of many instruments (Ewing & others 1957, Ewing & Press 1954(a), (b)). In 
addition, the group velocities obtained from studies of the direct wave agreed 
with determinations based on multiple passages of these waves around the earth, 
although the wave train had been strongly altered by its conspicuous dispersion. 
The agreement of these determinations shows that effect of the instruments 
on mantle Rayleigh waves is small compared to that of the dispersion. 

The essentially constant G wave group velocities were determined by Sat6 
(1958), exclusively by the measurement of differences in phase which were obtained 
by Fourier analysis of the arrivals G3, G1, and G4, G2. In this manner, the effects 
of the instrument and of the spectrum of the source are completely removed. 

The group velocities of long-period surface-wave signals provide a means of 
studying the mantle which is independent of body wave studies, and it is therefore 
quite important to get reliable surface wave group velocity data, as the derived 
mantle structures are sensitive to changes in the observed dispersion curves. 

After an investigation of the effect of seismographs on theoretical signals, 
calculations were made of the seismograms which would be produced by various 
instruments upon reception of the Gr signal from the New Guinea earthquake of 
1938 February 1, which had been used by Saté for his group velocity determination. 

Our calculations have shown that the instruments significantly influence the 
shape of pulse-like signals (see Figures 1 and 3) in cases of practical interest. Such 
changes in wave shape will cause measurable changes in the observed group 
velocity for signals with slight dispersion unless instrumental effects are removed. 
The modification of the G1 signal by the four pendulum-galvanometer systems 
described in Figure 4, which is similar to the distortion of the test signals in Figure 
3, indicates that the effect of the recording system on these pulse-like signals is 


* See Coulomb (1952), Di Filippo (1950), Ewing & Press (1956), Imamura (1937), Sat (1958), 
Sieberg & Gutenberg (1924), Wilson (1940). 
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quite important. Since this is the same G wave signal which, after elimination of 
the effects of the recording system and of the spectrum of the source, did not 
produce a maximum in the group velocity (Sat6é 1958), the earlier determinations 
of an essentially flat Love wave group velocity curve for long periods up to several 
hundred seconds are confirmed. This Love wave group velocity curve is useful 
for testing structures proposed for the upper mantle. An extensive programme of 
study of the upper mantle, using surface waves and body waves is nearing com- 
pletion at the Lamont Geological Observatory, and the results will soon be pub- 
lished (Landisman, Sat6 & Ewing). 

The present investigation of seismograph distortion of pulse-like ground 
motions indicates that, for best results, dispersion curves for pulse-like signals 
should be constructed from the seismograms of instruments whose magnification 
changes very gradually with period. In addition, this study emphasizes the fact 
that the transient response of instruments whose magnification changes rapidly 
with period can be important in studies of wave propagation. 


2. Method 


The transient response of various types of seismographs to pulse-like signals 
will be considered. The equations of motion of these systems, with suitable initial 
conditions, provide a complete description of their behaviour. 

One type of system is the linear strain as recorded by a galvanometer by means 
of a velocity pickup. The response, R, of a directly recorded strain rod, for har- 
monic components much longer than the rod, has been shown by Benioff (1935, 
1955) to be 


(1) 


for transverse waves, and 


(2) 


for longitudinal waves, 


where R = displacement of the free end of the horizontal strain rod relative to 
the earth beneath the free end, 
z = horizontal particle displacement, 
L = length of strain rod, 
c = horizontal phase velocity, 


a = angle between strain rod and the azimuthal direction of approach 
of the waves. 


Figure 1 in the study of Benioff & Press (1958) shows that in practical cases, 
changes of phase velocity of the order of 1km/s produce negligible changes of 
magnification. When phase velocity is treated as a constant and the time derivative 
of R, as obtained from a velocity pickup, is recorded on a galvanometer, the magni- 
fication is similar to that of a simple pendulum with the same period and damping 
as the galvanometer. The pendulum obeys the well-known equation 


¥+ 2egt + = — (3) 


where x = pendulum displacement with respect to the ground, 
z = absolute ground displacement, 


Oz 
Rr = — sina cosa — || 
c ct 
L oz 
= — — cos*¥a — 
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«9 = damping constant of the pendulum, 
9 = natural frequency of the pendulum. 


We may obtain the ground acceleration from a simple pendulum record (Bryan 
1936) by solving the equation 


= KjX+ Kgx, (4) 
where 


m= -I, kg = —2€9, and x3 = —ap?. 


The derivatives may be found numerically, using differences (Hildebrand 1956, 
Scarborough 1950). 

For a pendulum-galvanometer system, such as the Galitzin seismograph or 
similar instruments, where a mechanical pendulum drives a galvanometer through 
a velocity pickup, we may write the equations of motion (Coulomb & Grenet 
1935) in notation similar to that of Chakrabarty (1949): 


H+ + wo2x— oo) = —F = +, (5) 
6 + + — = 0, (6) 
where 
@ = galvanometer deflection, 
w = negative ground acceleration, 
€9 = pendulum damping constant, 
€y = galvanometer damping constant, 


«@o = pendulum natural frequency, 

(g = galvanometer natural frequency, 

oo = coupling of galvanometer to pendulum, 
6g = coupling of pendulum to galvanometer. 


It can easily be shown (Byerly 1942, Walter 1945) that the ground acceleration, 
= = —w, can be found if @, its first three derivatives, and first integral with 
respect to time are known. 


t 
#=-w= wel + + + | Ot (7) 
where t = time on record, 


Ko = 


Kl = —2(€9+ €g)/ 

K2 = + 4€0€9 + — 09, (8) 
K= €go")/ 0g, 

= 09. 


To obtain ground displacement we must solve 


t t 
| z= | ©) 
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When a linear strain record is to be compared with that produced by a pen- 
dulum-galvanometer system, it is best to use equation (4) to deduce earth accelera- 
tions and equations (5) and (6) to compute the comparison seismogram. ‘The 
equation (7) is only presented to show the method of intercomparison between 
various pendulum-galvanometer combinations. The computation of equation (4) 
will be referred to in future discussion as the inverse seismograph computation. 

To obtain a seismogram from the ground displacement or acceleration, one 
must solve the equations of motion (5) and (6). Many authors have contributed 
valuable studies concerning the effect of instrumental constants on the response 
of the instrument*, but it is often difficult to imagine what type of seismogram 
will correspond to a certain response curve and a given ground disturbance. 

We have adopted one of the Runge-Kutta methods, which has been derived 
by Gill (1951), to solve these differential equations numerically. We consider a 
system of n+1 simultaneous first-order differential equations 


d 
= = (10) 


where fo = 1, and yo = t. 


For such a set of equations, Gill has derived a numerical method of solution 
which is a fourth-order Runge-Kutta process, and in which, for each step of 
iteration over interval h, the error is of order h5. The process has the advantage that 
it may be started at any place in the iteration, and that the functions f; must be 
evaluated only at the ends of the interval and twice in the middle of the interval, 
which is convenient for interpolation or differentiation. 

The following set of equations, (11), which may be found in Gill’s study (his 


equation 26) are his solution for maximum accuracy. The equations are computed 
in the order shown, i.e. 


ki, then 7 for i = 1,...n, 


ya, then for i =1,...n, 


ya, then ga for =1,...n. 


The first of the two subscripts refers to the identification number of each member 
of the set of first-order differential equations, while the second refers to each of the 
four stages of the iteration process, which together constitute one complete step 
of interval h. When each of the sixteen equations, (11), have been computed for 
each member of the set of first-order differential equations, in the order mentioned 
above, the iteration has proceeded through one step of interval h. The q;; and r;; 
have been introduced to keep the errors, e(7,;), at each stage of the iteration process 
to a minimum. To preserve the accuracy from one step, /, to the next, it is neces- 
sary to make the gi of one interval equal to the g,9 of the next, and to round gio 
and gi4 in opposite directions. At the start of the first interval of iteration gio 


may be set equal to zero. With this introduction we may now consider the set of 
equations (11). 


* See Chakrabarty (1949), Coulomb & Grenet (1935), Eaton (1957), Hagiwara (1958), Matsu- 
moto (1958), Sutton & Oliver (in press). 
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kw = ---), ra = tho—qote(ri), 
ya = Gi = Got 3711 — 
ka = Afi(you, yu, ---), riz = [1— — Qe) + (2), 
ye = Yatne, giz = Gut 
Ria = hfi(yo2, -.-), ris = [1+ -V/(3)](Ri2— + 
= Gis = 3ria—[1 + 
kg = hfi(yos,y13; -.-), ria = + e(r14), 
via = Yistria, Gia = Fist 37u— 


A complete discussion of the method may be found in the original study by Gill 
(1951). 

Equations (5) and (6) were adapted to this method as follows: 

Let 


t=y, O= yz, 6 = yg. 
Then, y1 = ye, 


d 
or = = Silyo, V3, V4, w), 


dyo 
or 


= [—2e0v2— + covat 


Also, 


dy2 


f= — + ooyatw |. 


Similarly, ys = y4, or 


dy3 


Ja = [—2egya— oye), 


dys 


| fa = |. (6B) 


In order to find the acceleration, # = —w, from the displacement 2, the differ- 
entiated forms of the Stirling and Bessel central difference interpolation formulae 


106 
: 
(5A) 
or 
(3) 
Finally, 
or 
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to fourth-order differences (Hildebrand 1956, Scarborough 1950) were used to 
calculate w = 


I 
=- re — 0°0833 3333 322 + 1°3333 3333 321 — 2°5020 


+ 1°3333 3333 32-1 — 0°0833 3333 32-2], (12) 
and 
m= 6666 723 + 0°8125 O22 — 0°7083 3333 321 
— 0°7083 3333 320 + 0°8125 — 01041 6666 72-2], 
where wp is at the start, and w, is in the middle of an interval of iteration. The factor 
h is the interval of time between the various values of displacement, z; 20 is the 
displacement at time t; 21, 22 and 23 are later values of displacement, while z-; 
and 2z_2 are earlier values. If it is desired to use an acceleration input, the inter- 
polation formula becomes: 
wo = —%o, (13) 
and 
w, = —[+0°0117 1875 023 — 0°0976 5625 + 0°5859 375021 
+ 0°5859 3750%9 — 0°0976 5625 + 0°0117 1875 
At the start of the iteration process the central difference formulae, (12) and 
(13) must be replaced by Gregory-Newton forward difference formulae to fourth- 
order differences. Thus at t = 0, th, h and 4h we have, for a displacement input: 


=e jal + 09166 6666 721 — 46666 6666 723 + 95022 — 8:6666 6666 721), (12A) 


w= jal +0°2916 6666 724 — 16666 6666 723 + 4°25022 — 4°6666 6666 721), 


+ 09166 6666 725 46666 6666 724 + 816666 6666 722 
+ 2°9166 6666 721). 
jal + 6666 725 16666 6666 724 + 4°25023 — 4-6666 6666 722 
+ 1°7916 6666 721). 


For an acceleration input, we must use: 
= 0, (13A) 
= —[—0°0390 625024 + 0°2187 50%3 — 0°5.468 750%2 + 10937 504) ], 


— 
— [—0°0390 625025 + 0°2187 5024 — 0°5468 75023 + 1°0937 5022 


+ 0°2734 375041). 

For values of t 2 2A the central difference formulae are utilized. 

These equations, (4) for the inverse seismograph computation, and (11), 
(5A), (5B), (6A), (6B), (12), (12A), (13) and (13A) for the Runge-Kutta seismo- 
graph computation were coded for solution by the IBM 650 electronic data pro- 
cessing machine. 

A computation was made in which sinusoidal ground displacement signals 
were used as an input to a seismograph having the constants 7) = 308, €9 = 4°5wo, 
Tg = 1008, €g = wg, G9¢g = 0, which is one of the instruments in use at the 
Lamont Geological Observatory (Sutton & Oliver 1959). It will be referred to 
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hereafter as the 30-100 instrument. The steady-state response was the same as 
that calculated by Sutton & Oliver using the Laplace transform method of Chakra- 
barty (1949). See Figure 1, circles on curve c. 

A transient test signal was then made: 


I 2n 
t < t < t9+ 10s. 
2 20 


27 
cos—(t—10); te+10s < t < 308 
40 


I 2 
to+ 308 < t < t9+ 40s. 
2 20 


This ground displacement signal was used in a calculation with the same 30-100 
seismograph, and Fourier analysis of the input and output was performed (Matsu- 
moto 1958), with a Fourier analysis programme built for another study. 

The amplitude response to this transient 

GAT) = ——, where T = period, 
aT) (15) 
was then calculated, and is also shown in Figure 1 as curve c. The Fourier ampli- 
tude spectrum, A;(7), of the test signal S)(t) is also shown in Figure 1, as curve a, 


1-0 


Test signal S,(t) 
as ground displacement 


qT 


| To = 95 Tz = 90 


To Og = 0-3 


(Sign reversed) 


= 30 Ty = 100 
fo = 4°5Wo Wy 
Go Og 


qT 


(Sign reversed) 


Lj 


1 34 5 


time in minutes > 100 
T— 
Fic. 1.—Left: Seismograms produced by pulse-like ground motion, S,(t), on a 
narrow-band and a broad-band seismograph. 
Right: (a) Fourier amplitude spectrum of test signal S,(¢). 
(b) Fourier amplitude spectrum of test signal S,(t). Amplitudes 
times 0-1. 
(c) Amplitude response vs. period of 30-100 seismograph. Circles 


show results of steady-state tests, curve is the result of transient 
test. 
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and Ao(7) is the amplitude spectrum of the resulting galvanometer deflection, not 
shown. As can be seen in Figure 1, the transient response of the system is identical 
to that for steady-state signals. 

This same test signal was used to calculate the response curves shown in Figure 
2, and the input and output signals for the various instruments are shown in Figure 
3. The sign has been reversed where shown to facilitate comparison. Seismographs 
are normally connected so that an upward deflection of the trace implies that the 
ground moved up, east or north. 

The inverse seismograph computation programme was also tested with this signal. 
The resulting response curve for the linear strain recorder is shown in Figure 2, 
and the input and output time functions are shown in Figure 3. The resulting 
acceleration was then used as the acceleration input of the seismograph computa- 
tion. The resulting seismograms are shown in Figure 3. ‘The non-symmetrical 
ground acceleration shown in Figure 3 is that required by a 70 second pendulum 
in order to produce a seismogram equal to S;(t) (see equation (4)). By Fourier 
analysis of the 100-90 seismograms of Figure 3, and of the corresponding accelera- 
tion, the response curves for the 100-go instruments, shown in Figure 2, were 
obtained. 

In addition, Laplace transform calculation of the steady-state response at 
T = 60s was performed for the strain seismograph, and for the g5—go and 100~go 
instruments with og, = 0. These responses agreed with those found by Fourier 
analysis of the test signals, and are shown in Figure 2. 

From an assumed seismogram, S;(t), the acceleration was computed with 
equation (4). Then the computed seismogram which resulted, when a simple 
pendulum having the same constants received this acceleration as an input, was 
used to test the entire system of computation. The computed pendulum displace- 
ment signal was the same as the original input. Upon comparison, the first four 
leading digits of both signals were found to be identical. A similar test in which 
the test signal was first used as an acceleration, also reproduced the same accelera- 
tion to a corresponding number of digits. This means that the computed seismo- 
gram is the same as that which would be produced by the equivalent seismograph. 
The “noise” in the computed record is at least an order of magnitude less than the 
smallest deflections which can be plotted or observed. 

Another test signal, 


I 2m 
S2(t) = *|1—cos—| (to < t < to+20s) (16) 


was used as a displacement input to the 30-100 seismograph and the 95-90 instru- 
ment with coupling. The resulting signals are shown in Figure 1, and the Fourier 
amplitude spectrum of S2(t) appears as curve 6 in Figure 1. For signals which 
begin gradually, such as these, central difference formulae which start at tg # 0 
give equivalent results to forward difference formulae which start at to =o, 
to the precision required. Both methods were used and the results were equivalent. 

When it was clear that both computation programmes were working properly, 
calculations were made using an earthquake signal. A linear strain seismogram 
of the G1 signal from the New Guinea shock of 1938 February 1, as recorded at 
Pasadena, which was sampled every six seconds by Saté (1958) for his study of the 
G wave, was reduced to ground acceleration, and this acceleration was then used as 
an input for various approximations to the Gilman instrument discussed by 


— 100— 90 0, 0,~=0 


95— 9° = 0 


Eb — 100 — 90 Do Gy = 0:3(4£0tg) | fo = 
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Fic. 2.—Amplitude response vs. period for linear strain, T, = 70s, 


€, = w,, and four seismographs whose amplitude response changes 
rapidly with period. 
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Fic. 3.—Test input S,(t) used as ground displacement, left, and as 
strain seismogram, right. The 100-90 seismograms would be produced 
by the same acceleration that produces S,(t) on the linear strain. 
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Benioff & Press (1958). The results are shown in Figure 4. The overshoot dis- 
played by the transient signals in Figure 3 is evident in the resulting seismograms 
of Figure 4. The records in Figure 4 are those which would be obtained if all 
four pendulum—galvanometer instruments and the original strain seismograph were 
recording simultaneously. A comparison of the lower four seismograms shows 


Ground acceleration 
(Sign reversed) 


Linear strain 
T, = 70 sec 


time in minutes 
° 5 10 15 
i i Jj 


L 4 1 4 4 


Fic. 4.—Linear strain record, T, = 708, €, = w,, of the G1 signal 

from the New Guinea shock of 1938 February 1, as recorded at 

Pasadena (courtesy of Prof. H. Benioff) reduced to ground accelera- 

tion, above. The four seismograms produced by four approximations 

to the Gilman seismograph, below, display overshoot similar to that 
of the test pulses in Figure 3. 


that slight changes in the pendulum period have negligible effect, which is also 
shown by the response curves in Figure 2. The effect of changing the coupling 
from agsg = © to 0-3 (4€0€g) produces slightly increased distortion, for instruments 
whose pendulum and galvanometer periods are as close as these, but it is 
clearly evident that all four instruments have conspicuously distorted the seismic 


signal. 
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Figure 4 should be compared with Figure 8 in the study of Benioff & Press 
(1958). The record shown there, produced by a Gilman seismograph, is virtually 
identical with the computed seismograms in Figure 4. It was concluded by 
Benioff and Press that the group velocity of Love waves displays a maximum with 
a decrease in group velocity for periods greater than about two minutes, since the 
long period disturbance of the trace in Figure 8 of their study continued from the 
arrival of the G wave until the arrival of the Rayleigh wave train. ‘The original 
strain seismogram, shown in Figure 4 of this paper, was the same one which was 
used by Saté (1958) for a study of the G wave, in which the response of the instru- 
ment and the spectrum of the source were subtracted out by taking differences 
G3-Gi1, and G4-G2. These computations produced an essentially flat group 
velocity curve for long periods, which is reproduced in Figure 5. We are thus led 


6 


1000 


Wave Length — Km —> 


Fic. 5.—Phase velocity and group velocity of the G wave obtained 
by using G3 and G1 phases of the New Guinea earthquake (after 
Satd). 


to the conclusion that the long period disturbance recorded by the Gilman instru- 
ment after the arrival of a G wave is a purely instrumental effect, and is not related 
to the physical properties of the earth. 


3- Conclusions 


For best results, dispersion curves for pulse-like signals should be constructed 
from the seismograms of instruments whose magnification changes very gradually 
with period. This study also emphasizes that the transient response of instru- 
ments, whose magnification changes rapidly with period, may be important in 
studies of wave propagation. 

In the long period range up to several hundred seconds, there is no evidence 
for appreciable dispersion of the G wave, and earlier observations of an essentially 
flat group velocity curve remain valid. 
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The Determination of Virgin Strata Temperatures from 
Observations in Deep Survey Boreholes 


L. R. Cooper and C. Jones 
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Summary 


When a survey borehole is drilled downwards from the surface of 
the Earth the temperature is disturbed by the circulating drilling fluid, 
but fluid and adjacent strata gradually attain the local virgin strata 
temperature when drilling is stopped. Temperature equilibrium is 
approached most quickly near the bottom of the hole where the time of 
disturbance has been the least. Observations of the fluid temperature 
at this position may be made over a period of a day during a normal 
week-end break in drilling. Cylindrical heat flow theory is used in 
this paper to show that, with a knowledge of the thermal properties of 
the fluid and surrounding strata, the virgin strata temperature at the 
depth reached can be calculated from such observations. New and 
inexpensive equipment using a thermistor temperature measuring 
probe has been developed for making the observations of temperature 
and it has been tested successfully down to depths of 1 060 m (3 480 ft) 
in Hopton Pool Borehole in the Cannock Chase area. The paper also 
describes the equipment and gives the virgin strata temperatures for 
the Hopton Pool Borehole. It also gives observations of the tempera- 
ture of the borehole fluid at various depths taken on one day several 
days after the borehole drilling had ceased and shows that these could 
not be regarded as reliable measures of the virgin strata temperature. 


Symbols used in the text 


Virgin strata temperature. 

Borehole fluid temperature at a time tf = o. 

Borehole fluid temperature at time f. 

Temperature indicated by the probe (in °K). 

Radius of the borehole. 

Time. 

Thermal conductivity, density and specific heat, respectively, of 
the strata. 

Density and specific heat of the fluid within the borehole. 

Thermal diffusivity of the strata ( = k/pc) 

Electrical resistance of the thermistor in the probe (in ohms). 
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A,B Thermistor characteristic constants. 


H The reciprocal of the thermal contact resistance per unit area 
between the fluid in the borehole and the surrounding strata. 


t= Kt/a®, = 2pe/picr, = 9%), = ja. 


The units used throughout the report are in both the metric and foot-pound- 
second systems. 


1. Introduction 


When a coal mine is to be opened to bring coal from the warm strata deeper 
than goo m (2950 ft) it is necessary to know, as early as possible, whether the strata 
temperature is high enough to create an uncomfortable working climate. The 
earliest opportunity for obtaining information occurs when boreholes are being 
drilled downwards from the Earth’s surface to locate the depths and nature of the 
coal seams. During the progress of the borehole there are usually short breaks in 
the drilling during which temperature measuring equipment can be lowered into 
the hole. This paper describes new equipment which has been developed for 
measuring the temperature of the fluid at any depth in a borehole and also shows 
how the measurements very near the bottom of the hole can be converted to the 
virgin strata temperature at that depth. 

During the course of the work here described, it came to the authors’ notice 
that a similar programme was in progress for the Dutch State Mines but as far 
as is known to the authors no information has been issued on the results of the 
work and it is, therefore, not yet possible to make a comparison between the two 
cases. 

If a long hole of small diameter is quickly drilled vertically into virgin strata 
and the hole is immediately filled with a fluid which is cooler than the strata there 
will be a flow of heat into the fluid. At the bottom of the hole the heat flows 
from the sides of the hole and also from the base of the hole. However, it is not 
possible to lower a temperature measuring device to this position because the 
drilling rods require some time for withdrawal, especially from the deeper holes, 
and during this time the solid matter in the drilling fluid settles to the bottom. It 
is therefore not possible to measure the fluid temperature at the bottom and 
readings have to be taken a short distance above the base. In this position the 
heat flows into a short length of a cylinder of fluid and if the measuring position is 
twelve or more hole diameters above the bottom of the hole it has been shown 
(Blackwell 1956) that the heat flow is entirely radial towards the centre of the 
hole. For the special case of a survey borehole any theory associated with such 
temperature measurements must apply for different sizes of hole, for the diameter 
is decreased from about 24 to 12cm (g to 5 in.) as the hole proceeds downwards, 
and it must also take into account the thermal properties of the surrounding strata 
which change considerably, particularly in the coal measures. 

Normally, the effect of the drilling fluid on the strata temperature will pene- 
trate some distance into the strata during the long drilling process and Bullard 
(1947) has shown that some time would be required after drilling ceased for the 
recovery from the temperature disturbance. However, near the bottom of a 
rapidly-drilled hole, around which the strata have experienced little or no tempera- 
ture change, Bullard has indicated that the fluid temperature may approach quite 
near to the virgin strata temperature (V.S.T.) within as little as a day from the 
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end of the drilling. ‘This paper confirms this finding and shows that observations 
at even shorter intervals may be corrected to give V.S.T. 


2. Theory 


The radial flow of heat from an infinite region surrounding a cylindrical bore- 
hole to the fluid within it has been discussed theoretically (Jaeger 1956) and as 
long as there is no heat flow along the axis of the cylinder then 


exp( — tu?) du 
— (x — hu?) + [u Yo(u) — (x — hu?) ¥i(u) 


in which J,(u) and Yy(u) are the Bessel functions of order N of the first and 
second kinds. 

When a cool fluid is introduced into the borehole and its temperature rise is ob- 
served over a period of time it is more useful to study the ratio n = (6—9)/(@e—9) 


as this ratio approaches unity with increasing time. 
Let 
6-6, 
n= =I- = 1—F(h,x,7). 

For the case of a survey borehole in which the fluid is mainly water in good 
thermal contact with the wall of the hole then h = 0. Values of F(o, «, 7) have 
been evaluated (Jaeger 1956) and a graph of the values is reproduced in Figure 1 
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Fic. 1.—The function F (o.« 7). 
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where 
F(h, 2,7) = 
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for three values of «, namely 0-5, 1-0 and 1-5. For a borehole of known diameter 
containing water in a medium of known or assumed thermal properties at a higher 
temperature than the water, observations of the water temperature, @, (initially 
6 = 6) can be measured at a series of times, t. It is then possible to evaluate 
F(o, «, rt) from Figure 1, and hence n, for each of the times. The V.S.T., 4», 
can therefore be obtained for each time as 6, = 09+(1/n)(@ —09) and the values 
should be constant whatever the time. It is important to note that at a time t = 0 
the strata surrounding the borehole must be at uniform temperature. This 
method of temperature correction is the same as that used in other heat flow 
problems (Bullard 1954, Bullard, Maxwell & Revelle 1956). 


3- Experimental results 


That the above theory could be used very near the bottom of a hole was first 
established in the laboratory using an approximately quarter-scale model of the 
bottom of a deep borehole. The theory was then used in the determination of the 


V.S.T. surrounding the Hopton Pool Borehole in Cannock down to a depth of 
about 1 060m (3 480 ft). 


3.1. Examination of the theory by means of a laboratory scale model experiment 

A large block of sandstone of thermal conductivity 5 x 10~3cal cm~!s~1 (deg C)“! 
(0-00033 B.t.u. ft-!s—1 (deg F)-1, and about 1 m (yard) cube, was left to stand in a 
room until it has acquired a uniform temperature, the same as the surrounding 
air temperature. That this was so was checked by the readings of several thermo- 
couples, some placed near and on the faces of the sandstone and others inserted 
into very narrow holes drilled to various depths into the rock. A larger vertical 
hole of 4cm (14in.) diameter had been drilled downwards and centrally in the 
top of the block to reach just below the centre of the block. As soon as the tem- 
peratures in the block were uniform and steady the hole was filled with chilled 
slurry and a thermocouple was placed in the slurry eight hole diameters above the 
bottom of the hole. Figure 2 (full line) shows the temperature of the slurry over 
a period of time as it gained heat from the block. Table 1 shows the observations 
and the way they were corrected as described above to give the V.S.T. (chain 
dotted line in Figure 2) within 0-05°C of the measured V.S.T. which was 20°25°C 
(68-45°F). 

Satisfactory results were obtained after about half an hour of temperature 
readings. On a full scale borehole one would therefore expect to obtain similarly 
satisfactory results after about eight hours, the time changing with the square of 
the hole diameter (r = Kt/a?), It thus appeared quite possible to observe tem- 
peratures very near the bottom of a survey borehole for Jess than a day and so 
determine the V.S.T. Fortunately, it is possible to do this at the weekend break 
in drilling which usually lasts for about 24 hours. 


3.2. Field trials at Hopton Pool Borehole 


An opportunity to make observations of temperature in a survey borehole 
occurred in the Cannock Chase area where the National Coal Board is carrying 
out a programme of coal seam surveying. This borehole, known as the Hopton 
Pool Borehole and some 5km (3 miles) north-east of Stafford, commenced in 
1956 April and terminated at 1 060m (3 480 ft) at the end of 1957 September. 
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For about the first 150m (s5ooft), when the surrounding strata were water- 
bearing sand and pebbles, the hole was some 25cm (10in.) in diameter and was 
punched out. From 150 to 340m (500 to 1 100ft) the strata were red marl with 
subsidiary sandstone and the diameter of the hole was decreased to 18cm (7 in.). 
Between 340 and 450m (1100 and 1 5ooft) the strata changed to a mudstone 
which was impervious to water. Below 450m (1 500ft) the strata changed to the 
grey coal measures which were impervious to water and below 520 m (1 700 ft) they 
were predominantly mudstone with 21 coal seams between 520 and 860m (1700 
and 2840 ft), varying in thickness from 15cm (6in.) to 22m (7 ft 4in.). 


FC 
7°) Measured V.S.T. = 2025°C (68.45°F) Corrected temperature 


Observed temperature 


(Time hr) 
Fic. 2.—Laboratory-scale model experiment. 


Roller drilling began at about 150m (s5ooft) and core drilling commenced 
some 300m (1000 ft) deeper and continued to the bottom of the hole with a dia- 
meter just under 12cm (5in.). The hole was lined with a steel casing but as this 
never came nearer to the bottom of the hole than 30m (1ooft), and was often 
much further away than this, the probe used for temperature measurements was 
always surrounded only by strata. The fluid used for core drilling was a mixture 
of bentonite and water to give a specific gravity in the region 1-1 to 1-2, the higher 
figure being appropriate to the greater depths. Zero time was taken as the time 
when the drill rods were passing the depth at which observations of temperature 
were made. A short time elapsed while the rods drilled the short distance below 
the measuring point and some hours were required for all the rods to be removed 
from the hole and for the temperature probe to be lowered to its position near the 
bottom. This delay was timed as it entered into the correction to be made to the 
observed temperature. 

Before the field trials were commenced a theoretical estimate was made of the 
rise in temperature with time of fluid standing in boreholes of different diameters 
and surrounded by strata of different conductivities. As mentioned above, the 
Hopton Pool Borehole changed in diameter from 25cm (10in.) down to about 
12cm (5 in.) and the strata varied from sandstone to coal, giving a change in the 
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diffusivity in the ratio of about five to one. Figure 3 shows the ratio of the rise 
in temperature of the fluid to the initial temperature difference between strata 
and fluid obtained from calculations based on the theory given earlier. For a 
hole of 12cm (4}in.) in rock this ratio is greater than 0-98 at the end of 24 hours 
but some 160 hours are needed to achieve the same value for coal. The value of 
0-98 is selected here because for all values of the ratio greater than this the fluid 
temperature, 6, is within 0-5 °C (1 °F) of the V.S.T., 4, as long as 4) —8o is less 
than 28°C (50°F). 


(42in.) hole in rock 
24cm (9; in.) hole in rock 


Ring 
12cm (4: in.) hole in coal 


= 
Time (hr) 


Fic. 3.—Rise in temperature of the drilling fluid for 
different boreholes. 


Figure 3 also shows the effect of a larger hole in rock and, again, a longer time 
of some 100 hours is needed to achieve the ratio of 0-98. When the hole was large 
in diameter in the Hopton Pool Borehole at least 27 hours and up to 260 hours 
were available for observations and the difference between the virgin strata and 


drilling fluid temperatures was then the lowest, both being factors that reduced 
the temperature correction. 


3.2.1. The apparatus used for measuring the temperature of the borehole fluid 

The complete equipment is made up of the temperature probe, by which the 
borehole fluid temperature is measured, the measuring circuit by which the 
resistance of the thermistor in the probe can be determined, the cable connecting 
the two, the motorized winch and associated equipment for raising and Jowering 
the probe and cable and, finally, a depth measuring device on the cable. 
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Maximum and minimum thermometers were discarded in this project as it is 
necessary to observe changes in temperature over short times when it is neither 
desirable nor is there sufficient time to lower and raise the thermometers a suffi- 
cient number of times. Thermocouples were rejected as they would require 
compensating leads due to the steady change of temperature down a borehole, 
and the spurious voltages which occur in boreholes, which are referred to later, 
would upset the readings of the couples. Also taking into account abrasion 
resistance and weight of cable the use of thermocouples is not a practical proposi- 
tion. The choice of a temperature measuring device therefore rests between a 
thermistor and a wire resistance thermometer such as platinum, copper or nickel. 
Because the resistance of the latter is usually of the same order as the connecting 
cable it is necessary to use three or four leads to make full compensation for 
changes in cable resistance due to temperature changes. With a high resistance 
thermistor any change of cable resistance with temperature is negligible compared 
with the change in thermistor resistance and a cable with two conductors is all that 
is required. For the final design of temperature probe a type M.1512/300 ther- 
mistor made by Standard Telephones and Cables was used. Separate laboratory 
tests over a long period* justified the use of this type of thermistor. Figure 4 
shows a simplified section of the probe, of 2cm (?in.) diameter, the thermistor, 
having a copper base, being placed in direct thermal contact with the stainless 
steel case of the probe. The cable with its cable gland is detachable for transit 
purposes. The overall length of the assembled probe is 160cm (63 in.) and it is 
weighted with a lead filling to give a total weight of 2}kg (54 lb). 

The thermistor head was calibrated before and on several occasions during use 
against a standard thermometer in a bath of silicone oil, the thermistor resistance, 
R, and the absolute temperature, 7, being related by the law R = A exp(B/T) 
where A and B are constants. During the field work described in this report no 
significant change occurred in the thermistor calibration. 

The cable had to satisfy two separate requirements, one with regard to mechani- 
cal strength and the other with regard to its electrical properties. With a thermistor 
of about 100000 ohms at normal temperature a change of 0-05 °C (0-1 °F) in the 
temperature changed its resistance by about 250 ohms. ‘The series resistance of 
the cable wire, therefore, was satisfactory if it never changed by more than 100 
ohms. Similarly, a shunt resistance greater than 100 megohms changed the effec- 
tive measured resistance by less than 100 ohms. The cable had to withstand a 
safe load of 45 kg (100]b) (the maximum pull required to remove the probe from 
the mud) over and above its own weight at the full depth of the hole. The one 
finally selected was a special winch cable made by Johnson and Phillips and 
designated Type SS1A. It was made up of a central core of 7 tinned copper wires 
insulated with polythene and two armoured wrappings of galvanized steel wire. 
Its weight was 34kg/km (23 lb/1 ooo ft), loop resistance 128 ohms/km (39 ohms/ 
1000 ft), shunt resistance greater than 5 500 megohms/km (5 000 megohms/1 000 
yards) and its minimum breaking load was 410kg (goo lb). 

The measuring circuit consisted of a Wheatstone bridge (Figure 5), the probe 
thermistor forming one arm of the bridge, with a light-spot type suspended coil 
galvanometer having a sensitivity of 250mm/uA. When the bridge was in use it 
was found that stray voltages were observed of sufficient magnitude to deflect 
the galvanometer off scale. It is believed that these were caused by electro- 
chemical potentials developed along the length of the armouring which formed 

* made by L. R. Cooper. 
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Fic. 4.—Simplified section of probe. 
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one of the two conductors to the thermistor, and which was immersed in the 
corrosive and chemically active water in the borehole. By means of an additional 
battery and variable resistance an opposing voltage of the same magnitude was 
injected into the bridge circuit. Figure 5 shows this battery at Be with the added 

The external potentials were backed off by adjusting the 
potentiometer Rg to return the galvanometer spot to the mid point of its scale. 
During this operation the backing off battery key, Ks, and the galvanometer 
key, Ko, were closed, the bridge battery key, Ki, was open. The final balance was 
then obtained, with Ke and Ks still closed, by adjusting the decade resistances R; 
until there was no change of deflection of the galvanometer on opening and closing 
K;. The resistance of the thermistor probe was then read off from the dials of Rj. 


Galvanometer 


Fic. 5.—Bridge circuit. 
Ri decade 10 to 211 100 


That the added circuit did not affect the accuracy can be seen by noting that 
the bridge balance was obtained by observing the effect on the circuit only of 
battery B;. The other e.m.f.’s played no part in this operation. Considering the 
circuit with all e.m.f.’s other than B; replaced by short circuits the condition for 
balance is that of a four-arm Wheatstone bridge. It is then seen that the resistance 
Rg is practically equal to the resistance of the circuit R5,¢,7,8, at all values of Rg and 
thus the main circuit is unaffected. The overall accuracy of the bridge in measur- 
ing temperature was better than 0-03 °C (0-05 °F). 

A builder’s power winch was specially adapted for the field work by the manu- 
facturers. It consisted of a two-stroke engine which drove, via a four-speed gear 
box (three forward speeds for raising the probe and one reverse for lowering), a 
pair of cork covered rollers. The drum which carried the cable normally rested 
against brake pads but, on depressing the control handle, the drum moved off the 
brakes on to the cork driving rollers. There was thus a “‘dead man’s handle” 
effect for quick stopping. The weight of the winch with cable was about 136kg 
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(300lb). The pull was 160kg (3501b) at a minimum speed of 8 m/min (25 ft/min). 
A depth indicator was added to show the amount of cable paid out, its smallest 


indicated unit being 7}.cm (3 in.). Figure 6 (Plate I) shows the apparatus assembled 
for use. 


3.2.2. The borehole temperature measurements 


Because of the probable interchange of water between the strata and the bore- 
hole down to about 150m (sooft) no figures are included to this depth in this 
paper although the apparatus was given preliminary tests at the depth of 50m 
(160 ft). The first reliable figures, capable of correction in accordance with the 
theory, were obtained below the water-bearing strata. 

It was not possible to measure the temperature of the drilling fluid as the drill 
passed the depth at which the probe was subsequently positioned and it was 
therefore not possible to give a value to %. However, it is not necessary to know 
this value as long as several observations of the probe temperature are measured 
over a period of a few hours. It was shown in Section 2 of the paper that 


= —F(h, 2, 7)] 


and so it is possible to plot the observed values of @ against the corresponding 
values of {1 —F(h,x,r)}. By least squares the best straight line through the 
observed points can be obtained and this will have a gradient and an intercept from 
which the virgin strata temperature can be evaluated. This technique, when 
applied to the results of the last hour of the scale model experiments in Table 1, 
gave the value of 20°3°C (68-5°F) for 6) compared to the measured value of 
20°25°C (68-45°F). The technique was adopted with the borehole temperature 
measurements. 

At the depth of 156m (51o0ft) the probe gave a steady temperature of 10°7°C 
(51°3°F) for some 16 hours of observations. As these observations were com- 
menced ten days after the drill rods had been withdrawn the correction to be made 
to the observed reading was negligibly small (Figure 3) and the value of 10-7°C 
(51°3 °F) was accepted for 4). 

At the depth of 302m (ggoft) the probe was lowered near to the bottom of 
the hole ten hours after the drill rods had drilled past this depth. ‘The observed 
temperature was steady at 13-3 °C (55-9°F) for 22 hours and the analysis of the 
observations gave 6, as 13.3 °C (§5-9°F) there being no correction on this occa- 
sion. 

When the probe was lowered on a later occasion to 500 m (1 640 ft) the bottom 
of the hole was smaller in diameter and was passing through grey coal measures. 
Figure 7, curve (a), gives the observed temperatures for a period of two hours and 
the correction required to give the V.S.T. of 17-45 °C (63-4°F). Curve (b) gives 
corresponding figures for the borehole when at a depth of 786m (2 580ft) and 
6, = 29°05 °C (84°3°F) and curve (c) gives 0, = 31°6°C (88-9°F) at 878m 
(288oft). In these three cases the probe readings were commenced some 11 or 12 
hours after the drill rods had drilled past the depths at which observations were 
made. The hole was just under 12cm (5in.) in diameter and the strata were 
predominantly mudstone with coal layers. 

The last opportunity for taking temperatures in the Hopton Pool Borehole 
occurred when the hole ceased at 1060m (3 480ft). Before the probe could be 
lowered to the bottom it was necessary for Schlumberger to carry out two traverses 
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Fic. 6.—Apparatus assembled for use. 
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to the bottom of the hole and for this the borehole had to be flushed with freshly 
made slurry. Unfortunately for the probe measurements the hole was con- 
tinuously flushed, without drilling, for just over two days prior to the withdrawal 
of the rods. Some cooling of the strata must therefore have occurred and it was 
expected that the results would not be so satisfactory on this occasion. Observa- 
tions were therefore carried out for as long as was possible, some 17 hours, and the 
maximum observed reading was 37°45 °C (99°4°F), the rate of increase of tem- 
perature at the end of the observations being 0-5 °C (1°F) in 10 hours. It was 


V.S.T. 316°C (889°F) 


(c) 878m (2880 feet) 


VS.T. 29.05°C (84:3°F) 


(b) 786 m (2580 feet) 


Temperature 


V.S.T. 17.45°C (63.4°F) 
(a) 500 m (1640 teet) 
Observed 


° 2 
Time (hours) 


we 


Fic. 7.—Examples of the Hopton Pool Borehole temperatures. 


considered that on this occasion the theory could not be applied. It does not give 
a temperature change of the order observed, which was unusually high, and the 
preliminary flushing of the borehole should have reduced the heat flow below the 
theoretical value. It is certain, however, that the maximum observed temperature 
was below V.S.T. and from previous experience at other depths it is likely that 
this temperature should be increased by at least 1°C (2°F) to give a value of 
39 °C (102°F), which is probably not far removed from the V.S.T. at 1060m 
(3 480ft). Such a value fits the general trend of the graph of temperature and depth 
for the coal measures (Figure 8) which one would not expect to change unless 
there was a marked change in the strata conductivity. 
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3.2.3. An underground measurement of virgin strata temperature at Brereton 
Colliery 

During the drilling of the Hopton Pool Borehole an opportunity was taken to 
make one measurement of V.S.T. in the coal face at Brereton Colliery using a 
different technique to be reported elsewhere (Harris & Jones). Brereton Colliery 
gives the nearest underground access for comparison with the borehole figures. 
It lies about 16km (10 miles) south-east of Stafford and a little more than this 
distance from the borehole. A V.S.T. of 13.4°C (56-1 °F) was observed at 317m 
(1 ogo ft). 


4- Discussion of results 


The results described above are summarized in Table 2 and Figure 8 shows the 
V.S.T. plotted graphically against depth of cover. 


Table 2 
Virgin strata temperature at Hopton Pool Borehole 


Depth ‘Temperature 
Maximum observed 


°F 
10°7 
13°3 
16°55 
27°4 
30°9 87-6 
37°45 99°4 


Brereton Colliery 


Also shown on Figure 8, as crosses, are the temperatures observed with the 
probe at different depths as it was raised from the bottom of Hopton Pool Bore- 
hole on the last visit. The probe was rested at each depth for some minutes until 
its reading was constant. It will be seen from Figure 8 that the temperatures 
obtained during this vertical traverse through the borehole fluid lie below the 
V.S.T. line at greater depths and above it towards the surface. One would expect 
this sort of result for during the long time of drilling a borehole the strata towards 
the bottom are cooled by the drilling fluid and those at the top are therefore 
heated, a phenomenon generally appreciated by users of temperature logging 
equipment when used in seam survey (Leonardon 1936). Such readings are there- 
fore of only small value when accurate V.S.T. is required unless, of course, the 
hole can be left undisturbed for some time (Bullard 1947), two or three weeks if 
accuracy within 0-5 °C (1°F) is required. When times of a few hours only are 
available in which to determine V.S.T. it is therefore necessary to observe 


| 

m ft °C °F 
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temperature as quickly as possible very near the bottom of the hole, as described 
in this paper, repeating the observation as the bottom reaches various depths. 

Figure 8 also indicates that the geothermic gradient is not constant and varies, for 
the Hopton Pool Borehole, from about 50 to 26 m/deg C (go to 48 ft/deg F) between 


Sand and [Red marl 
| pebbles FF Mudstone | Coal measures 


100 5 


O Bottom of the bore hole 
+ Vertical traverse 


Brereton Colliery 


800 1000 m 


1000 1500 2000 2500 3.000 3500 feet 
Depth of cover 


Fic. 8.—Virgin strata temperatures at Hopton Pool Borehole. 


the depths of 300 and 1060m (1000 and 3 sooft). It is of interest to compare 
this change with the calculated value when the thermal conductivity of the strata 
changes, the Earth’s radial heat flow remaining the same. Table 3 gives the geo- 
thermic gradient assuming the Earth’s heat flow (Benfield 1939) is 10~* cal cm~*s~1 
(3-7 x 10-6 B.t.u./ft?s) and using estimated values of the thermal conductivity of 
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the strata between the depths tabulated. The conductivities are as realistic as 
possible, short of actual measurement in situ, and are based on a knowledge of the 
geology of the strata obtained during boring and on laboratory measurements of 
strata thermal conductivities (Benfield 1939, Godridge 1956, Lomax 1957). 


Table 3 
Change in geothermic gradient due to change in conductivity of the strata 


Depth Strata Strata conductivity Calculated 
temperature geothermic 
gradient 


°F cal/emsdegC B.t.u./ft hr deg F m/deg C ft/deg F 


13°35 56-0 
0-006 1°45 73 133 
14°0 57°2 
0-005 5° 
16°45 61-6 
00033 0-80 33 
18-0 64°4 


0°003 3° 
36-05 96-9 


The change in the gradient is very similar to that observed, confirming that 
the gradient is considerably affected by the lower conductivity of the strata in 
the coal bearing measures. 


5. Concluding remarks 


(a) By observing the temperature of the drilling fluid very near the bottom of 
a borehole for no more than a day it is possible to calculate the virgin strata tem- 


perature of the surrounding strata. A method of calculation has been developed 
for this purpose. 


(b) The thermistor temperature probe, the resistance bridge, the motorized 
winch and ancillary equipment, all newly designed and inexpensive, are satisfac- 
tory for making the observations, having been tested on many occasions quite 
successfully during the drilling of a borehole to a depth of 1060m (3 480ft) at 
Hopton Pool in the Cannock Chase area. 


(c) The temperature of the borehole fluid at various depths, taken on one 
day several days after the borehole had finished drilling, is below virgin strata 
temperature towards the bottom of the hole and above it at the top of the hole. 
Borehole fluid temperature can become a direct measure of the virgin strata tem- 
perature only if the hole has been left standing for some time. 


(d) At the Hopton Pool Borehole the virgin strata temperature changed with 
depth in a manner which was not linear, from 10°5°C (51°F) at 156m (510ft) 
to about 39°C (102°F) at 1060m (3 480ft). The geothermic gradient changed 
from about 55 to 26m/°C (100 to 48ft/°F) for the same change in depth, due 


probably to the reduction in the thermal conductivity of the strata in the lower 
coal measures. 


302 99° 
351 1150 
473 1550 

525 1720 
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Rayleigh-Wave Dissipation Functions in Low-Loss Media 


J. Ross Macdonald 
(Received 1959 March 3) 


Summary 


A simple formula is derived relating the specific dissipation 
factors for Rayleigh, compressional, and shear waves in relatively 
low-loss materials. Calculations and a table are presented which 
allow an unknown specific dissipation factor to be obtained directly 
from the two known factors for the other types of waves. Results 


apply for any realizable ratio of any two of the three elastic phase 
velocities. 


1. Introduction 


The surface Rayleigh waves which travel in a solid half-space in addition to 
body compressional and shear waves are of particular importance for earthquake 
detection at large distances, where they predominate over three-dimensional 
body waves. As Press & Healy (1957) have pointed out, it is of interest to relate 
the absorption coefficients of Rayleigh waves to those of shear and compressional 
waves for materials showing energy dissipation and dispersion. ‘These authors 
have derived a useful expression for this purpose for low-loss media which we 
shall here simplify, slightly generalize, and reduce to numerical applicability. 

It will prove most convenient to deal with the specific dissipation factor, 1/Q, 
for each type of wave, where 27/Q is the ratio of energy dissipated per stress 
cycle to the peak energy stored (Knopoff & MacDonald 1958). The quantity 


1/Q is a useful quantity for comparison with experiment because it is independent 
of the detailed mechanism of energy dissipation. 


2. Dispersion and absorption 


Whenever dissipation occurs, the resulting dispersion makes it necessary to 
distinguish between three kinds of wave velocities: the elastic velocities in the 
absence of dissipation, the corresponding phase velocities with dissipation, and 
the complex velocities. Using the subscripts A, B, and C to denote quantities 


relating to compressional, shear, and Rayleigh waves, the elastic velocities may be 
obtained from 


va = [(A+2p)/p}', (1) 
vp = (2) 
= 16[1 — (3) 
where p is the density and A and yw are Lamé’s moduli. 
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A damped harmonic wave travelling in the x direction involves the factor 
exp[iwt —yx], where y=«+i8 and « and 8 are the attenuation and phase factors 
respectively of the wave. For a given type of wave, the phase velocity V is given 
by w/B8 while the complex velocity is V’ = w/(—#y) = w/k. Further, when dis- 
persion is present, the real equation (3) is replaced by the corresponding equation 
involving the complex velocities V;’. Following Press & Healy, this complex equa- 
tion can be used to determine a relation between 1/Q,4, 1/Op, and 1/Q, in low- 
loss materials. It is assumed in the following that these three quantities refer to 
waves all having the same period. 

From their definitions, it can be shown that 


Vi/Vii = (4) 

1/Qi = — (5) 

Now if the dissipation is relatively low so that 1/O <o-1, (5) can be expanded and 
simplified to yield, with negligible error, 

Vi/ Vi,’ | —i(1/2Q;), (6) 


where terms in (1/2Q;)? have been neglected. To the same order of approxima- 
tion, it can be shown that for either the Voigt dissipation model considered by 
Press & Healy or for loss arising from creep (Lomnitz 1957) V4~va and Vg~ vp. 
Taking Vce~vc as well and defining 


-a=(vc/va)*, b=(vc/vz)*, (7) 
the complex-velocity Rayleigh equation becomes 


1—i(1/2Qc) 1—i(1/2Qc) 1—i(1/2Qc) 


When (8) is expanded retaining terms up to 1/Q; only, one obtains an equation 
from the real part identical with (3) and an equation from the imaginary part 
involving the 1/Q;’s. By using (3) to simplify the latter, or by taking imaginary 
parts of the logarithms of (8), one finally finds 


I I I 
= (1-m) + m—, (9) 


Oz On 


a(2—b)(1~5) 
a(2 —b)(1 —b) —b(1 —a)(2 — 36) (10) 


Although it is obtained more generally, Equation (9) is a simplified but equivalent 
form of the equation obtained by Press & Healy (1957). 

The result (g) is a simple linear equation connecting the 1/Q;’s and is valid 
to high approximation when none of them exceed 0-1. This equation can be used 
to obtain any unknown 1/Q; given the other two and knowledge of m. It is 
interesting to note that 1/O is a more sensitive measure of loss than is dispersion 
of the velocities which, to the present approximation, is negligible. 

In order to make (g) useful, the quantity m, which is a complicated function 
of v4, vg, and vc, has been calculated, using a digital computer, as a function of 
vp/va and is presented in Table 1 and Figure 1. Table 1 also gives vc¢/vg as a 
function of vg/va, and the resulting curve is shown in Figure 1 with its ordinate 
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on the right. Now, if any two of the elastic phase velocities are known, the third 
can be obtained using the Figure or the Table. Then, after the corresponding value 


of m is obtained, the unknown 1/Q; can be obtained using (9) and known values 
of the other 1/Q,’s. 


Table 1 


Numerical results 


0°95531 ° 
0°95516 3°2411 X10-4 
0-95469 1-3233 
095389 3°0823 
0-95271 5°7570 x 10-* 
O°95112 9°6002 x 10-3 
©°94903 1°5008 x 1072 
0°94633 2°2597 X 10-2 
0°94286 3°3347 10-2 
0-93837 4°8851 x 10-2 
093253 7°1789 x 10-2 
0-92476 0-10685 
O-91419 0-16257 
0°89937 0°25512 
0°87785 0°41586 
0°70455 


T 


T 


0 86 
°7 


Fic. 1. The quantities m (log scale) and vc/vg versus vz/va. 
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It will be noted that the abscissa extends only from o to 0-7. The ratio vg/va 
may be expressed for isotropic material as 


(vs/va) = (11) 
where v is Poisson’s ratio. For an incompressible body v = 4 and v4 = o. 
The minimum value of v is zero; thus the maximum value of vg/v,4 is 2+. The 
quantity m increases continuously beyond this physical limit, however, and reaches 
g°1 for vg/v4 = 0°95 and infinity for vg/v4 = 1. 

When one or more of the 1/Q;’s are appreciably greater than o-1, the exact 
complex equivalent of (3) must be solved using (4) and (5) in their exact forms. 
Solution of the resulting cubic equation with complex coefficients will be diffi- 
cult but is usually unnecessary since the condition 1/Q;< 0-1 is usually well met 
in most situations of physical interest. 

Finally, it should be pointed out that for vg/v4 > 0-567008 there are three real 
roots of (3) rather than the single one which appears for smaller vg/v4. As 
Jeffreys (1952) has shown, however, it is only the smallest of the three roots which 
is of physical significance for Rayleigh wave transmission. 
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Research Note 


An Earth Model and the Reduction of Gravity 
Observations for use in Stokes’s integral 


J. C. Bhattacharji 
(Received 1959 February 23) 


Summary 


Application of Dr Hunter’s Earth model in the reduction of 
gravity observations for use in the Stokes integral is discussed. An 
attempt is made to make it clear that there are practical difficulties in 
utilizing the new surface anomaly designated by Dr Hunter as g—y x, 


in the Stokes integral, for the purpose of determining ultimately the 
figure of the Earth. 


Dr Hunter derives his Earth model from the actual Earth by smoothing its 
topography without change of total mass. What is actually required is to compute 
the weighted average values H of the actual heights A of the land topography (or 
—o°615 times the depths of the ocean bed) for a sufficient number of control 
points, each average height being over a region of 150km radius. ‘These computed 
values enable us to draw the contours of average height, which actually define the 
Earth model. It is thus “bounded by a surface of gentle and slowly varying 
slopes” never exceeding 1/100 even in the worst cases of rugged topography, 
which have but negligible effects on gravity. It is, therefore, considered advan- 
tageous to reduce an observed gravity value to the surface of the Earth model, 
that is, to a point having its height equal to the weighted average height over a 
region of 150km radius, by means of the usual Bouguer reduction and then to 
compare it against the theoretical value of gravity reduced to the same height by 
the free-air method in order to obtain the gravity anomaly designated by Dr Hunter 
as g—Yz- The disturbing effect of the more irregular part of the Earth’s actual 
topography, having been removed as above, the surface anomaly g—y, becomes 
necessarily more regular than either g—y, or g—yg and, therefore, as repre- 
sentative of the unseen anomalies of masses as any other isostatic anomaly derived 
by heavy reductions. Moreover, the surface of the Earth model being very 
smooth compared with that of the real Earth, the gravity values once reduced to 
the spherical surface of the layer of thickness H, limited only to 150km radius, 
do not require appreciable corrections to be applied for any departure of this 
spherical surface from the actual surface of the Earth model, while by reducing 
them to the smoothed surface of the latter the anomaly g—y, becomes comparable 
with similar anomalies anywhere in the world without requiring any assumptions 
of internal densities of the Earth. This is surely a great advantage. So far there 
is no difficulty. 
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But things are entirely different with the geoid. The geoid is affected signifi- 
cantly by the superficial transfers of masses within and beyond the limit of 
150km radius, and is, therefore, not a level surface of the Earth model. ‘‘Corres- 
ponding to the Earth model will be a co-geoid—a level surface of the same 
potential as that of the geoid.” Now according to the Helmert’s second method of 
condensation, all matter external to the co-geoid can be condensed into a surface 
lying in the same co-geoid by reducing the gravity value on the surface of the 
Earth model by the free-air method using the normal rate of change of gravity 
with height above the co-geoid, along with the Gelandereduktion added to it, which 
in the case of the Earth model is apt to be inappreciable. The vertical separation 
between the actual co-geoid and the ideal co-geoid after condensation is also 
negligible. So after the transfer of the external matter, either co-geoid encloses 
the attracting matter of the Earth model, thereby reducing the co-geoid to a 
bounding equipotential surface, and the corresponding gravity anomaly on being 
reduced to the same surface, may be used correctly in Stokes’s integral to give the 
height of the co-geoid above the reference spheroid. It is thus essential to compute 
the vertical separation between the geoid and the co-geoid (or the consequent 
indirect effect of gravity) which is required for reducing not only the gravity anoma- 
lies to the co-geoid after condensation, but also the height of the co-geoid derived 
directly from Stokes’s integral to the required height of the geoid above the 
reference spheroid, by considering only the topographical effect due to the super- 
ficial transfers of masses. Now the mass-transfers for the Earth model being only 
superficial, visible features having been removed to fill visible hollows, the shift of 
the centre of mass and of the rotation axis are likely to be negligible but the deforma- 
tion of the geoid, on the other hand, is quite significant and what is still worse, 
the actual computation of the corresponding vertical separation between the geoid 
and the co-geoid due only to the topographical effect of the mass-transfers, is 
very much complicated in the absence of considerations of the effect of isostatic 
compensation, which in fact considerably reduces the large topographic effect 
increasing normally with distance. It is also quite clear, from Lambert and 
Darling’s tables for computing the deformation of the geoid on the Hayford 
system based on the hypothesis of complete isostatic compensation, that the 
resultant geoidal effect of distant masses or of their displacements is an easily 
computable small quantity inappreciably affected by any possible error of estima- 
tion of distant topography or of their displacements and is therefore obtainable 
with sufficient accuracy, though the effects of topography and compensation 
considered separately are large quantities significantly and almost equally affected 
by a slight change in position or amount of distant masses and are, as such, bound 
to be extremely inaccurate, even if computed with all necessary care and with 
difficulties involved in such reductions. This is surely an unfortunate position but 
cannot be avoided if our object is to utilize Dr Hunter’s surface anomalies in 
determining ultimately the height of the geoid above the reference spheroid. This 
will be made more clear from the following example of the topographic and the 
isostatic reductions of the deformation of the geoid for a Himalayan station about 
8-5km in height, on the basis of the hypothesis of complete compensation at a 
depth of 113-7km. 

In this example the mean heights of zones refer to the estimated mean 
heights of Hayford zones in the actual Earth. These are likely to differ from the 
corresponding values on the Earth model, the discrepancies becoming more pre- 
dominent in regions of rugged topography but decreasing slowly with distance. 
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To take account of the topographic effect of such discrepancies is again not only 
difficult but also practically unobtainable with the desired accuracy even after heavy 
reductions. 

Thus it would seem clear that the promise held by this new surface anomaly 
£—Yx in the field of geodesy has got its limitations and as such may not lead us 
far in the realization of our objective of determining the figure of the Earth. 


10 Raipun Road, 
Dehra Dun (U.P.), 
India. 


1959 February 18. 
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= effect 
in cm 
7°00 94 77 
6-60 82 65 
6-00 165 124 
5°80 227 159 
6-02 408 255 
5°70 619 325 
4°35 1466 546 
3°73 1700 360 
3°20 2464 259 
18 2°86 677 43 
17 2°79 802 40 
16 2°76 988 38 
15 2°69 1240 36 
14 2°64 1633 32 
13 2°54 40°79 43 
12 2°44 4494 27 
11 2°32 6113 I 
10 0°450 
3-036 6317 12 
9 —0°075 
4691 12 
8 4727 27 
7 —1-812 
0-915 741 5 
6 —2°452 
0°487 577° 27 
5 —3°475 
0-240 16811 74 
4 a —26282 100 
3 —4°019 
48996 | 200 
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I —2-807 —8005 38 
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Letter to the Editors 


Magnetic “ Cleaning ”’ 
E. Irving 


On the basis of demagnetization experiments with 14 specimens from one rock 
formation in southern France, As & Zidjerveld (1958) conclude (p. 318) that 
“from these demagnetization experiments it follows that the pole determinations 
that are based on measurements of remanent magnetization in rock from which 
the disturbing magnetization of the present Earth’s field has not been removed 
should be treated with reserve.” 

Now the mean direction (D = 216, ] = —24) of the specimens studied by As 
and Zidjerveld after removal of the secondary component (that produced by the 
present Earth’s field) is in broad agreement with those published earlier by Roche 
(1957) who gives D = 210, J] = —16 as the mean of five pyromeride specimens 
and D = 185, J = —13 for three dolerite specimens collected from the same rock 
formation. No reference is made to this earlier work. Moreover the pole position 
they give (47N 144E) agrees very well with that obtained from work on other 
Permian formations of Europe without magnetic “cleaning”: Exeter traps 43 N 
164E (Creer 1955); Montcenis sandstone 38 N 162 E, St. Wendel sandstone 45 N 
175 W, Nideck rhyolites 43 N 168E (Nairn 1957); Oslo porphyries 46N 165 E, 
(Rutten & others 1957); Lavas 36 N 175 E, sediments 37 N 163E and intrusives 
32 N 174E of Ayrshire (Du Bois 1957, Armstrong 1957); Upper Permian of N.W. 
Russia 45 N 178E (Khramov 1958). The agreement is very good when it is 
remembered that each pole determination is subject to errors of the order of 10°, 
and that each of the formations is not exactly contemporaneous and so small polar 
movements between the times represented by each are not excluded. These new 
results given by As and Zidjerveld do not cast suspicion on what has been done 
previously; on the contrary, they confirm it. In fact this agreement between 
many determinations from Permian rocks in Europe would seem to be a most con- 
vincing confirmation of procedures in palaeomagnetism. 

The techniques of magnetic “‘cleaning’”’ of rocks are providing the means of 
dissecting out the unwanted secondary components which often mask the primary 
components in the natural state. Up till now only a few per cent of rocks have been 
of use, but now our range of vision may be increased many times. The mood is 
one of optimism for the future, not reserve about the past. 


Dept. of Geophysics, 
Australian National University, 
Canberra, Australia. 


1959 April 3. 


References 


Armstrong, D., 1957. Nature, 180, 1277. 
As, J. A. & Zijderveld, J. D. A., 1958. Geophys. J., 1, 308. 
140 


| 
| 
| 


Letters to the Editor 


Creer, K. M., 1955. Thesis, Cambridge. 

Du Bois, P. M., 1957. Adv. in Phys., 6, 177. 

Nairn, A. E. M., 1957. Bull. Soc. Geol. France, 7, 721. 

Khramov, A. M., 1958. Palaeomagnetism and Stratigraphic correlation, 187. 

Roche, A., 1957. C.R. Acad. Sci. Paris, 244, 2952. 

Rutten, M. G., Everdingen, R. O. & Zijderveld, J. D. A., 1957. Geol. Mijnbouw, 
19, 193. 


Note on Progress in Geophysics 


Dynamic Processes of Abyssal Sedimentation: Erosion, 
Transportation, and Redeposition on the Deep-sea floor 


Bruce C. Heezen 


Summary 


The concept of a quiet, nearly motionless, and certainly currentless 
deep-sea floor lacking the processes of erosion so prevalent in shallower 
seas was well ingrained in scientific thought until discoveries of the past 
decade revised this serene and static picture. The new concept is a 
more dynamic one in which submarine landslides, raging turbidity 
currents, internal tides, deep-sea current scour, and submarine post- 
depositional solution and alteration play a part in shaping the sea floor, 
in disturbing and nurturing its fauna and which produce sedimentary 
structures and distributions formerly thought to be the exclusive mark 
of shallow water deposits. ‘The concepts involved are still in a critical 
state of flux with many conflicting viewpoints as to the processes 
involved and their relative importance. The development of the major 
concepts is traced and the principal variant viewpoints are briefly 
reviewed. 


Part 1. Early ideas 


The Naked Abyss. ‘“The geological clock may strike new periods; its hands 
may point to era after era; but, so long as the ocean remains in its basin—so 
long as its bottom is covered with blue water—so long must the deep furrows and 
strong contrasts in the solid crust below stand out bold, ragged, and grand. 
Nothing can fill up the hollows there; no agent now at work, that we know of, 
can descend into its depths, and level off the floors of the sea.” (Maury 1855.) 

Thus M. F. Maury eloquently summed up the view generally held a century 
ago. Remembering the then newly discovered deep-sea oozes he added: 

“But it now seems that we forgot the oceans of animalculae . . . whose remains 
in vast multitudes sink down and settle upon the bottom.” 


The Unremitting Snowfall. The concept of the great unremitting snowfall 
which modified and then replaced the original picture of a naked abyss was stated 
clearly by Rachel Carson: 


“When I think of the floor of the deep sea, the single overwhelming fact that 
possesses my imagination is the accumulation of sediments. I see the steady, 
unremitting, downward drift of material from above, flake upon flake, layer upon 
layer—a drift that has continued for hundreds of millions of years, that will go 
on as long as there are seas and continents.” 
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Large elements of truth are found in both of these seemingly contradictory 
concepts which is why both have so long endured. To be sure, sedimentation 
is so slow over vast areas of the ocean floor that the original ‘‘grand and imposing” 
relief mentioned by Maury is, over broad areas of the sea floor, still to be observed 
only slightly modified by subsequent deposition (Heezen, Tharp & Ewing 1959). 
Furthermore, the success which students of sediments have had in reading that 
“epic poem of earth history” contained in the successive layers of sediment 
further attests to the general validity of the concept of the unremitting snowfall 
(Phleger 1952, Ericson et al., in press, Broecker, Turekian & Heezen 1958). 
However, neither of these concepts tells the complete story; nature is considerably 
more complex than either model. 


A Revolutionary Discovery. John Murray announced in an obscure footnote 
added to a discussion on deep-sea sediments published in 1912 that: 


“From his examination of the samples collected on board the Gauss, 


Philippi believed that stratification on the sea floor of today is not the exception 
but the rule.” 


This revolutionary discovery led to hypotheses of uplifted and deeply down- 
faulted, mobile blocks to account for each successive thin stratification. The most 
perplexing of the stratifications were those which involved an alternation of deep- 
sea ooze or clay on one hand, with sand and silt of a type previously found only 
near the shore. As long as only-a few such deposits were known (Andree 1920), 
the problem was not faced straightforwardly, but when, after World War II, 
dozens upon dozens of examples were brought up in deep-sea piston cores, the 
problem reached tremendous proportions. 


The Atlantic Cable. Maury (1855) discussing the practicability of the pro- 
posed telegraph cables wrote with confidence: 


“It is an established fact that there is no running water at the bottom of the 
deep sea. The agents which disturb the equilibrium of the sea, giving violence 
to its waves and force to its currents, all reside near or above its surface; none of 
them have their home in its depths.” And, therefore, he concluded that if a deep- 
sea cable were laid it ‘‘would lie in cold obstruction, without anything to fret, 
chafe, or wear, save alone the tooth of time”. That this clear statement also con- 
tains a large element of truth is shown by the great success of the telegraph cables 
which have broken infrequently under blue water. However, not long after the 
first cable was laid, the rude fact that this picture did not tell the whole story 
became very apparent. 


A Rude Awakening. In an early report Benest (1899) wrote: 


“Unmistakable evidence, offered in some cable repairs in recent years, has 
disclosed novel causes of rupture, and has extended our knowledge of the hidden 
dangers that beset the telegraph wire after it has been assigned to its oozy bed.” 


Early investigations of the distribution of earthquake epicentres were made 
with the practical aim of avoiding earthquake belts in laying out future cable 
routes, for it was early shown that earthquakes and cable breakage were closely 
related. In his classic paper on ‘“‘Suboceanic Change” Milne (1897) documented 
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many cases where a positive relationship was proved. Benest, writing at about the 
same time, noted the seasonal breakage of cables along certain continental slopes, 
and thus did not strongly favour the seismic cause of cable interruption. He 
noticed that many earthquakes failed to produce cable failures; this, together with 
the seasonal occurrence of breaks led him to suspect a connection with the seasonal 
floods of nearby streams. His appeal to submarine underground river outlets, a 
misguided effort, was revived more recently by Douglas Johnson (1939). 

The mode of emplacement of deep-sea sands and the cause of deep-sea cable 
breakage remained independent problems for nearly three-quarters of a century. 
It has been only in recent time that these two phenomena were related to a 
common cause. 


Part 2. Turbidity currents 


The Unifying Concept. ‘Turbidity currents, the concept which provided a 
unified explanation for these previously unrelated problems was not a brand new 
theory; but, like so many other concepts in science, it was resurrected after being 
ignored or rejected for over fifty years. Interest in turbidity currents originally 
stemmed, not from their relationship to the origin of deep-sea sands, nor from their 
responsibility for the breakage of cables, but from the explanation of another related 
mar on the picture of a motionless pristine abyss. 

The discovery of submarine canyons had been made originally by wire sound- 
ings; but these giant gorges carved into the submerged flanks of the continental 
blocks became a real enigma when, as the result of careful surveys made with the 
then newly-developed echo-sounder, it became clear that canyons could be found 
in large numbers wherever surveys were made on the continental margins. As 
an explanation, the first and most popular theory involved the drowning of sub- 
aerially eroded valleys (Spencer 1903, Shepard 1948). However, this concept 
rapidly lost favour when it was demonstrated that submarine canyons reached the 
very bottom of the ocean basins instead of being limited to the upper parts of the 
continental slope. R. A. Daly (1936) resurrected the original concept of F. A. 
Forel who had proposed in 1885 that the sublacustrine canyon of the Rhéne was 
cut by currents of sediment-laden water flowing in along the floor of Lake Leman 
(Lake of Geneva). Inspired by Daly’s eloquent appeal Kuenen (1937) and Bell 
(1942) successfully produced turbidity currents in tanks. Grover & Howard 
(1938) actually observed weak turbidity flows in Lake Meade Reservoir (Gould 
1951). 


Muddy Water. ‘Turbidity currents produced by the early experiments flowed 
at relatively slow speeds, did not erode, and were, in fact, little more than halting 
underflows of muddy water (Kuenen 1937). Stetson & Smith (1937) believed that 
the density contrast between turbidity currents and the superadjacent water was 
so small that a marine turbidity flow could leave the bottom and flow as a turbid 
interflow when it reached a depth where the density of the water exceeded that of 
the turbidity current. In later experiments the use of higher density mixtures 
yielded greater velocities (Kuenen 1950). Although high density turbidity cur- 
rents had been effectively demonstrated in Kuenen’s experiments, scepticism and 
lack of confidence in the scale factors prevented a general acceptance of the process 
when projected to full scale. The concept of low-density turbidity currents which 
leave the bottom to form interflows was also adopted by Menard & Ludwick (1951). 
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Although they concluded that submarine landslides ‘‘are the only likely sources 
of dense turbidity currents thus far suggested”’ they had serious reservations, and 
felt that such turbidity currents resulting from slides ‘‘may be tenuous rather than 
dense”. They also concluded that “turbidity currents resulting from suspended 


sediment in rivers must be rare’. Concerning the generation of turbidity currents 
along the shore they concluded: 


“If they exist, marine turbidity currents produced by wave agitations, or any 
other stirring agent, probably do not have high effective density.” . . . 


‘“‘No means exist for holding the current in place until the suspended sediment 
becomes very concentrated; therefore, it is unlikely that very dense turbidity 
currents can be generated even by violent wave and current agitation.” . . . 


‘*. .. turbidity currents may flow on the ocean surface or within the ocean due 
to marine density stratification. . . . Surface flows and interflows may be important 
geologically because they can spread sediment over wide areas.” 


For this reason they conclude that “marine turbidity current deposits may be 
difficult to identify”’. 


The concept thus projected was not appealing and most oceanographers and 
geologists still favoured the picture so eloquently phrased by Maury (1855). 

The case for modern high density, high velocity turbidity currents was finally 
proved by a study involving the explanation of both the breakage of submarine 
cables and the origin of deep-sea sands. 


Full-scale Experiment. In 1952 Heezen and Ewing concluded that the: 


“*. .. events associated with the Grand Banks earthquake of 1929 November 18 
may be considered as a full-scale experiment in erosion, transportation, and deposi- 
tion of marine sediments by a turbidity current in which the submarine telegraph 
cables served to measure its progress and to give tangible evidence of its force’’. 


Following this earthquake which shook the continental slope south of New- 
foundland, all the submarine telegraph cables lying downslope (south) of the epi- 
central area were broken in sequence from north to south. According to Heezen 
and Ewing’s explanation each successive cable was broken by a turbidity current 
which had originated as a series of slumps in the epicentral area. The current 
travelled down and across the continental slope, continental rise, and ocean basin 
floor and continuing far out onto the abyssal plain, reached a limit well over 450 
miles from its source area on the continental slope. 

On the continental slope where the bottom gradient is 1:10 to 1:30 the velocity 
of the current exceeded 50 knots; and out on the abyssal plain where the gradient 
is less than 1:1500 its velocity, although diminished, still exceeded 12 knots. 
Subsequent exploration of the area with the piston corer revealed a one-metre 
thick, uppermost bed of graded silt containing shallow-water microfossils (Heezen, 
Ericson & Ewing 1954) overlaying abyssal lutite. The thickness of the bed was 
nearly identical with that calculated by Kuenen (1952) on the basis of the slopes 
and observed velocities. 

This bold challenge to the concept of the motionless, pristine abyss was imme- 
diately attacked. Shepard (1954) based his objections largely on the claim that 
the velocities of 10-50 knots were larger than could be expected, but offered no 
more serious objections. Kullenberg (1954) misled by some erroneous soundings 
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claimed that Heezen and Ewing’s explanation required the turbidity current to 
run uphill and over a range of hills, an ability which, he correctly pointed out, 
could not be attributed to turbidity currents. However, subsequent exploration 
confirmed Heezen and Ewing’s belief that the alleged mountain range was the 
result of poor contouring of bad soundings. In fact, a survey by Heezen, Ewing & 
Ericson (1954) revealed nothing but a smooth abyssal plain in the indicated posi- 
tion of Kullenberg’s mountains. 

One of the strongest critics of the turbidity current concept is Professor Hans 
Pettersson (1954), the well-known leader of the Albatross Expedition. 'To explain 
the occurrence on the equatorial Atlantic abyssal plains of graded, deep-sea sands 
which contained “‘twigs, nuts, and the bark of dicotyledonous trees” he offered 
the following alternate to turbidity current deposition: 


“In the event that a large island, harbouring vegetation, and with a fairly exten- 
sive shelf covered the Mid-Atlantic Ridge northwest of St. Paul’s Rocks and 
became submerged during a catastrophe of seismic-volcanic character a few hun- 
dred thousand years ago, material like that found in the Albatross cores might have 
become distributed over the adjacent sea bottom.” 


Abyssal Plains. Following World War II continuously recording echo- 
sounders were first extensively used in deep water. One of the first discoveries was 
that, while most of the deep-sea floor is covered by hills and mountains, large areas 
in the deepest portions of the ocean are covered by abyssal plains (Heezen, Ewing 
& Ericson 1954). The origin of these vast flat areas posed a real problem until the 
concept that turbidity currents were filling the oceanic depressions offered a ready 
solution (Heezen, Ewing & Ericson 1951). Pettersson and some of his followers 
including Koczy remained unconvinced, and to this day prefer to explain the 
abyssal plains as the result of great outpourings of lava (Koczy & Burri 1958). 

During the exploration of the abyssal plains, long narrow steep-sided channels 
were discovered which cut into the abyssal plains. These features, called mid- 
ocean canyons (Ewing & others 1953), or deep-sea channels (Dietz 1953, Menard 
1955) were indeed an enigma to those who believed that the plains are the result 
of the outpourings of lava. Koczy & Burri (1958), in an attempt to explain the 
occurrence of features normally ascribed to erosion and deposition in an area 
where these authors deny the existence of such processes, offered the suggestion 
that mid-ocean canyons were, in fact, fault-bounded rift valleys, and that their 
natural levees were instead the sides of tilted fault blocks. This explanation says 
nothing of the graded sand beds found in the canyon floor, the graded silt beds 
found on the levees, nor of any of the other streamlike characteristics of the North- 
west Atlantic Mid-Ocean Canyon (Ewing & others 1953). Recent magnetic pro- 
files across this Mid-Ocean Canyon failed to reveal an anomaly, a fact which can 
be added to the list of objections to Koczy and Burri’s theory. 


Submarine Natural Levees. Submarine natural levees have been observed 
along the continental slope canyons of southern California (Buffington 1952), as 
well as along the mid-ocean canyons of the Atlantic (Ewing & others 1953), the 
Indian Ocean (Dietz 1953), and the deep Pacific (Menard 1955). The largest deep- 
sea natural levees observed to date are those of the Congo Submarine Canyon 
(Heezen & others 1957) which in a depth of 2 250 fathoms are 100 fathoms high and 
over 20 miles wide. It is safe to conclude that the submarine natural levees are 
of nearly the same origin as the subaerial ones. When turbidity currents overflow 
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the tops of the canyon walls they suddenly deposit a considerable portion of their 
load as they abruptly slow down and disperse. 

Since the natural levees of the Northwest Atlantic Mid-Ocean Canyon are 
nearly identical in form with those bounding the deeper parts of the Congo Sub- 
marine Canyon, it seems unlikely that they are the result of a radically different 
process as proposed by Koczy & Burri (1958). 


Deep-sea Sands. Ericson, Ewing & Heezen (1952) presented a preliminary 
account of their study of turbidity-current sediments deposited in the abyssal 
plains of the North Atlantic. Some of the most dramatic evidence came from two 
cores taken from the narrow trench abyssal plain which lies at about 4 580 fathoms 
at the bottom of the Puerto Rico Trench. These cores consisted of abyssal red 
clay at their base; but each contained a graded bed of calcareous sand, silt and clay 
which included fragments of the algae Halimeda. Since living Halimeda requires 
sunlight this detritus must have originated near the shore. 

Over 200 cores have been obtained from the abyssal plains of the North 
Atlantic. Each core taken in an abyssal plain contained either sands or silts, many 
with shallow water Foraminifera. Most cores from the flanks of the Mid-Atlantic 
Ridge, the abyssal hills and the oceanic rises contained only pelagic deposits, 
although a few at the bases of islands or from the crest zone of the ridge contained 
coarser material (Ericson & others, in press, Heezen, Tharp & Ewing 1959). 


Incredible Coincidence. After the initial shock wore off, a general acceptance 
of deep-sea turbidity currents as an important process was found in most quarters. 
Pettersson (1954) admits the possibility of their existence but concludes that: 


“The present evidence hardly affords clear proof that turbidity currents have 
played any dominant part in sculpturing the ocean floor. Both their extension 
over the deep ocean basin and their velocities of steep slopes appear to have been 
over-estimated by Heezen and Ewing.” 


Shepard (1955) now feels that turbidity currents transport sediments, but 
denies that they have the ability to erode. Unable to accept the high velocities 
indicated by the Grand Banks sequence of cable brakes, Shepard (1954) suggested 
that an unrelated series of local and independent submarine landslides had broken 
the cables. He wrote: 


“Is it not possible that slides developed more slowly on the gentle outer slopes 
and built up pressure gradually against the cables until they broke?” “‘. . . other- 
wise the sequence could be explained as at least in part a coincidence and in part 
due to the slower development of slides at an increasing distance from the epi- 
centre.” 


Jones (1954) also believed the breaks were: 


“’. . much more likely to have been caused by slumping than by the thin 
stream of material carried in a turbidity current (flowing at) an incredible velocity”. 


Spontaneous Liquefaction. Encouraged by the objections to the turbidity cur- 
rent hypothesis raised by Shepard and others, but not satisfied with their alternate 
explanations, Terzaghi (1957) offered a new explanation for the cable breakage 
which followed the 1929 Grand Banks Earthquake. He supposed that the sedi- 
ment in the epicentral area and over the whole area of cable breakage to the south 
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existed in a metastable state. He supposed that the earthquake shock caused spon- 
taneous liquefaction over the epicentral area and that the apparent velocity of the 
turbidity current was instead: 


**... the velocity of the advance of the boundary of the liquefied portion of the 
sediments underlying the ocean floor and not the velocity at which the slide material 
moved over the floor’’. 


The only fact he attempts to explain is the velocity indicated by the sequence 
of cable breaks. He is neither concerned with the problem of how the sediment 
came to be in the supposed metastable state, nor why the sequence of breaks was 
downslope and not in all directions from the epicentre. Neither is he concerned 
with the emplacement of the shallow water sediments in deep water nor the grading 
of the deposited beds. It is not readily apparent why the passage of such a wave of 
liquefaction should break submarine cables with at least 10 per cent slack. One 
might think that the cable would be required to sink thousands of feet into the 
sediment before breakage would occur. Although it cannot be denied that the 
mechanism he suggests may locally affect sediments, the weight of evidence seems 
to definitely point away from Terzaghi’s theory of propagation of spontaneous 
liquefaction as an explanation of the Grand Banks cable breaks. 


Undetected Earthquakes. A series of cable failures on the floor of the western 
Mediterranean essentially similar in distribution and order to the Grand Banks 
sequence followed the 1954 Orleansville earthquake (Heezen & Ewing 1955). 
These events further supported the interpretation of the Grand Banks turbidity 
current and provided additional evidence of the mode of emplacement of the Medi- 
terranean deep-sea sands which Bourcart (1953) had previously interpreted as 
turbidity-current deposits. 

Bourcart & Glangeaud (1956, 1958), however, objected to the interpretation 
claiming that turbidity currents had broken only the cables on the steep slopes. 
To explain the later breaks out on the nearly level ocean floor they suggested that 
undetected or at least unlocated aftershocks had caused the breaks. Apparently 
impressed by Kullenberg’s objections to the computed velocities for the Grand 
Banks turbidity current, they sought an alternative interpretation not requiring 
such high velocities. Also misled by bad soundings, they denied the existence of 
an abyssal plain in the area. Precision soundings, however, have confirmed the 
existence of the Balearic Abyssal Plain. 


The Big Creep. Lombard (1956), referring to the transportation of ancient 
sediments, and Koczy & Burri (1958), referring to the present deep-sea environ- 
ment, favour as an alternative to turbidity currents, a slow creep which might 
require hundreds or thousands of years to traverse a distance that a turbidity cur- 
rent would cross in only a few hours. That slumps, slides and other slower and 
more massive slope failures occur in the ocean is clear from the evidence of 
topography, of sediments and of sedimentary structures. But the results of such 
slow movements are quite different from the graded, well-sorted deposits laid down 


by turbidity currents travelling at express train speeds over the ocean floor. 
Koczy (1954) notes that: 


“Near the shelf the surface is often undulating because of the presence of the 


sediment sliding down the continental slope. This movement can be assumed to be 
very slow.” 
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This slow creep, although certainly of importance over wide areas, cannot be 
considered as a substitute for turbidity currents, but rather as an additional and 
independent process (Kuenen 1956). 


‘Elevator Tectonics.”” A concept which might be described by the term “‘eleva- 
tor tectonics” (or perhaps “‘lift faulting”) has been often employed in the explana- 
tion of stratified deep-sea deposits (Cushman 1941, Landes 1952). The latest 
proponents of this concept are Greenman & Le Blanc (1956) who invoke a 10000 
foot subsidence of the entire central portion of the Gulf of Mexico at the close of 
the Wisconsin (11 000 years ago) to account for the emplacement of grey silts now 
found on the floor of the Sigsbee Abyssal Plain. These grey silts have been ex- 
plained by Ewing, Ericson & Heezen (1958) as the results of turbidity-current 
deposition of sediment contributed by the Pleistocene Mississippi River. Cores 
from low knolls rising above the abyssal plain contained only slow, uniformly 
deposited pelagic sediments. According to Greenman and LeBlanc’s explanation, 
the knoll core stations would have stood more than 150 metres above sea-level until 
commencement of Recent time. Radiocarbon measurements on the pelagic deposits 
from the knolls indicated post-glacial rates of deposition of 3cm per 1 000 years, 
while the Pleistocene rate in the same cores ranged from 10 to 20cm per 1 000 
years. In contrast, the sediments cored from the abyssal plain showed post-glacial 
rates of 10cm per 1 000 years and glacial rates of 60cm per 1 000 years. 


Archipelagic Aprons. Broad, smooth, deep-sea aprons encircle many oceanic 
islands. Ericson, Ewing & Heezen (1952) reported graded beds of coarse, shallow 
water detritus from the Bermuda Apron. More recently, Menard (1956) referring 
to similar features observed in the southwest Pacific, coined the term ‘‘archipelagic 


apron’. He reasons that while turbidity currents may have smoothed their surface, 
the aprons owe their existence to extensive outpourings of lava. Since their 
volume exceeds that of the pedestal they surround, he reasons that turbidity cur- 
rent deposition cannot account for their entire volume and consequently that a 
gentle rise toward the island pedestals must have existed before the smoothing 
process started. Such features lacking a smooth surface can be detected near some 
islands and seamounts but they are excluded by their roughness from the defini- 
tion of archipelagic aprons. 


Tsunami. Bucher (1940) formerly thought that tsunami waves set up by the 
Grand Banks Earthquake had caused the later sequence of cable breaks. Indeed 
there seems to be a relationship between tsunami and cable breaks; however, they 
are probably not cause and effect but both result form the motion of the slump and 
ensuing turbidity current. Heezen & Ewing (1952) concur with Gutenberg’s 
(1939) view that “principal cause of each tsunami is submarine slumping.”” How- 
ever, the turbidity currents themselves probably are also capable of generating 
tsunami. That tsunami could not be the direct result of earthquakes was proved 
by the fact that not all large submarine earthquakes produce tsunami. In order to 
prove that turbidity currents and slumps alone could produce a tsunami, we clearly 
need to find a tsunami and a related series of cable failures not preceded by an 
earthquake; or to disprove the concept, to discover a turbidity current lacking an 
associated tsunami. 


Hydrodynamic Considerations. ‘There was considerable feeling that the velocity 
and extent of the Grand Banks turbidity current were over-estimated. Robert 
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Stoneley (1957) investigating the problem from the viewpoint of hydrodynamics 
concluded that “‘the hypothesis would be much more workable if the current had 
a thickness of only a few metres rather than over 100 metres suggested by Kuenen’s 
calculations”. He found that the question of 


. .. whether the motion of the current will persist over distances of the order 
of 1000km without being damped out by viscous forces . . . hinges on the viscous 
drag of the stream at high Reynolds numbers. Theoretical and experimental 
work on arrested saline wedges suggests a much lower drag coefficient than is 
appropriate for slower speeds, and the extrapolation to high Reynolds numbers 
again suggests a thinner turbidity current of vertical thickness about 5 metres’’. 


The thickness resulting from Kuenen’s (1952) calculations had also seemed 
excessive from a geological point of view since turbidity-current deposits were not 
found on slight elevations above turbidity-current blanketed abyssal plains. 


Staublawinen. Careful study of snowslides may provide instructive analogies 
for the understanding of the behaviour of slumps and turbidity currents. Roher 
(1954) who fired smail charges in the air a half-metre above the snow surface, 
triggered snowslides in a small Alpine valley (Grat des Sazmartinhorns). He cast 
lines of coloured dye on the snow and by observing with both pressure gauges set 
in the snow and coloured motion picture photography, was able to differentiate 
between turbulent and laminar slides. He found that in laminar slides the velocity 
did not exceed 32km/hr. The sides of a laminar slide were sharp; and when it 
stopped, the slide left a cone of sediment on the valley floor. Under certain snow 
conditions the laminar flows, after moving about 100 metres, would undergo a 
quick transition (usually during less than 30 metres of movement) and flow in a 
turbulent manner. In the same valley where laminar flows travelled at 32 km/hr, 
turbulent flows of 114km/hr were observed. The velocity of turbulent flows timed 
in a variety of locations ranged from 85 km/hr to 305 km/hr. The deposits of dyed 
turbulent flows were spread out evenly on the valley floor with the dye evenly 
mixed through the deposit. The transition along the sides of turbulent flows was 
a gradual one. In one case the suction produced in the air by the passing of a fast 
moving flow pulled men from one side into the slide. In Staublawinen the snow 
mixes with air, while in turbidity current sediments mix with water; and although 


the two may not be directly analogous, their comparison may help in the better 
understanding of high velocity turbidity currents. 


Requirements of Turbidity Currents. Turbidity currents first of all require sedi- 
ment, just as rivers require water, without it they do not exist. Turbidity currents 
require a trigger mechanism; this can be supplied by earthquakes, hurricanes 
impinging on the shore, high bedload discharge of rivers; and lacking a trigger 
mechanism, turbidity currents probably can occur after long continued deposition 
simply by slope failure resulting from over-steepening of a depositional slope. 
Finally, turbidity currents require a slope. The optimum conditions for the genera- 
tion of a turbidity current probably would be a large body of fine, well-sorted sand 
triggered by an earthquake on a steep slope. 

Heezen (1956) studying cable failures from various parts of the world has 
found evidence for several different trigger mechanisms. Besides earthquakes, 
hurricanes have, in at least one instance, produced a turbidity flow off Cape 
Hatteras. One of the most interesting types of trigger mechanism is that related 
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to the high bedload transport at the mouths of rivers. Potentially, such a source 
should have a greater geologic importance since both the trigger effect and the 
source of sediments are supplied by the river. In such an area repeated flows from 
a single source could build up a submarine abyssal cone, and at the same time 
repeated currents down the same path on the continental slope could erode a canyon. 
Significant periodic changes in depth have been observed at the heads of sub- 
marine canyons by Shepard (1948), Zenkovich (1958), Heezen (1956) and others. 

It was found that cable failures in the submarine canyon of the Magdalena 
(Heezen 1956) occurred after prominent changes had taken place in the depth and 
shape of the river mouth. These changes in turn had occurred only during the 
highest flood stages of the river. A similar pattern was also observed off the Congo. 
Both rivers have a common characteristic in their lack of a modern subaerial delta. 
The delta sediments in both cases are to be found in the alluvial cones at the 
mouths of their submarine canyons; and indeed, exploration in both areas has 
disclosed graded beds containing a large proportion of terrestrial organic remains 
as well as tree leaves and wood. A dramatic fact associated with the 1935 cable 
failure off the Magdalena was recorded when grass, still green, was brought up 


wrapped around the damaged cable which lay in 750 fathoms, 12 miles seaward 
from the river mouth. 


Abyssal Life and Petroleum. The occurrence of carbonaceous deposits in the 
deep-sea alluvial cones led to speculation on the effect of this material on the 
nutrition of the abyssal fauna (Heezen, Ewing & Menzies 1955, 1957). It was 
reasoned that turbidity flows would have two radically contrasting effects on the 
abyssal fauna. Rapid deposition of a thick bed of sediment over the abyssal floor 
must have a tremendous destructive effect, smothering by burial much of the 
indigeneous fauna. In contrast, the large windfalls of organic debris brought to 
the deep-sea by turbidity currents must support a greater productivity than the 
normal, slow, pelagic deposition. 

Repeated turbidity currents from an organic-rich source build an organic-rich 
submarine deposit containing both the transported organic debris and the trapped 
and entombed abyssal fauna. 

Another important effect of turbidity currents on petroleum occurrence is 
independent of the bottom fauna. Orr, Emery & Grady (1958) found that while 
the content of chlorophyll derivatives in sediment immediately below turbidity 
current deposited sands was anomalously high, that of the sands themselves was 
anomalously low. The turbidity currents which brought in the sands were appar- 
ently derived from a source low in organic matter and thus in this case only the 
blanketing effect of the turbidity current deposit is seen. 

Turbidity currents flowing into an anaerobic basin (such as the Cariaco Trench, 
Heezen & others 1958), may provide sandy source beds and impermeable clay cap 
beds for future oil fields. In this case the source of the oil is the planktonic 
remains. It seems likely that much of the oil in the productive Maracaibo basin 
was formed and trapped in this manner. 

Since the geologist’s key to the past is the present, the implications of sub- 
marine turbidity currents as regards paleogeographic reconstructions are tre- 
mendous. No longer can scattered shallow water fossils be considered evidence of 
shallow-water deposition, nor does evidence of strong currents indicate shallow 
depths. Even ripple marks, long considered an absolute proof of shallow water 
origin, have been found as deep as 2 000 fathoms. 
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Modern Geosynclines. ‘Turbidity currents are gradually filling the modern 
deep-sea trenches. Turbidity-current smoothed trench (abyssal) plains have now 
been reported from every trench which has been sounded in detail. Trenches 
bound most of the circumference of the Pacific Ocean, but in the Atlantic such 
marginal topographic trenches are exceedingly rare. However, wherever seismic- 
refraction studies have been made at the base of the continental slope a deep trough 
of lower velocity material ascribed to unconsolidated sediments and sedimentary 
rock has been found paralleling the base of the continental slope. These former 
trenches were probably filled largely by turbidity current sediment. They consti- 
tute in all probability, the modern geosynclines which quietly await the process 
of orogeny which through deformation and uplift, may eventually turn them into 
folded mountains. 

The effect of turbidity currents over a long period of geologic time is to fill in 
the hollows of the deep-sea floor by building vast abyssal plains, and to erode 
canyons into the continental slopes. Of course, if the rate of subsidence of the 
marginal trenches exceeds the rate of deposition, broad basin floor abyssal plains 
would not form since all turbidity current deposits would be trapped in the 
trenches. The effect of turbidity currents in the deep sea produces a topography 
essentially similar to that produced by periodic downpours in desert mountains. 


The coalescing alluvial fans resemble the continental rise, and the playas resemble 
the abyssal plains. 


Continued Acceptance. Turbidity current deposits have been described from 
the modern basins of sedimentation off southern California by Gorsline & Emery 
(1959). These modern turbidity-current deposits have been compared by them to 
those described by Natland & Kuenen (1951) and others from the filled sedimen- 
tary basins of southern California. Considerable attention has been given to these 
deposits since the analogous deposits in the filled basins are productive of petroleum. 

Menard (1955) has employed the turbidity current hypothesis to explain the 
smooth topography, and the canyons, and the deep-sea channels of the North 
Pacific. Locher (1952, 1954) and Phleger (1951) who studied the lithology and 
fauna of the deep-sea sands from the equatorial Atlantic became convinced of 
their turbidity-current origin. Bourcart (1953) and Dupliax & Calleux (1956) 
described deep-sea sands from the abyssal plains of the Mediterranean which they 
identified as the deposits of turbidity currents. The concept has found an in- 
creasingly greater application in palaeogeographical reconstruction and a host of 
recent papers have appeared which we will not attempt to review here. It will 
suffice to say that Kuenen, his associates, and his students as well as many other 
independent workers have had considerable success in identifying ancient turbidity- 
current deposits throughout the world. 


Over-enthusiasm. Unfortunately, the turbidity-current concept has been 
occasionally overworked even in regard to recent deposits. By invoking turbidity- 
current action Rigby & Burckle (1958) attempted to explain the emplacement of 
thin beds (millimetres thick) of fresh-water diatoms found in cores of mid- 
equatorial Atlantic eupelagic sediment. ‘This explanation is no less ridiculous 
than the theory proposed by Malaise (1957) who believes them to be deposits laid 
down in ponds during a recent and complete emergence of the Mid-Atlantic 
Ridge. The Congo River at present carries fresh-water diatoms hundreds of 
miles offshore. Kolbe (1957a, 1957b) points out that they are also blown hundreds 
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of miles to sea by the wind. To emplace diatoms on the Mid-Atlantic Ridge as 
suggested by Rigby & Burckle (1958) would require turbidity currents to climb 
over 3000 metres in 500 miles travelling all the way over an intensely rugged 
bottom, an ability not attributable to turbidity currents. Rigby and Burckle 
do not doubt that surface currents, or perhaps the wind, can transport diatoms as 
far as the core stations on the Mid-Atlantic Ridge, but they do not believe, for 
some unspecified reason, that ocean currents could deposit and concentrate beds of 
fresh-water diatoms. 

The geological record contained in the ancient sediments is sufficiently complex 
that other misapplications will undoubtedly occur. 


Part 3. Erosion at depth 


A New Surprise: Abyssal Ripples. When photographs of the deep sea floor 
were first taken in appreciable numbers the frequent discovery of ripple marks and 
scour marks was simply astounding. In 1948 a well-known physical oceanographer, 
upon seeing a picture of ripple marks at 400 fathoms was so amazed that he 
seriously suggested to the writer that the ripples were fossil and had been formed 
before a local subsidence of the sea floor. Ripples were found in many localities 
inaccessible to downslope movements. It is improbable that turbidity currents 
are related in any general way to the processes forming the ripples seen in bottom 
photographs. Actually, the ripples should not have come as a surprise since early 
studies of sediments taken from submarine peaks and in deep straits in depths 
ranging down to 2000 fathoms had revealed winnowing of the finer particles 
(Murray & Renard 1891). In fact, Agassiz (1888) early recognized that current 


scour could be effective to depths down to at least 650 fathoms. He concluded 
that: 


**.. . the bottom of the Gulf Stream along the Blake Plateau is swept clean of 
slime and ooze, and is nearly barren of animal life’. 


Weber, in 1900, later noted coarse sand and rock bottom to depths of 1 500 
metres in the passages between the Indonesian Islands. But concerning these 
findings Grabau (1913) concluded that: 


“It is probable that we have here a prevention of sedimentation by the removal 
of the particles before they reach the bottom, rather than any effect of eroding 
work of the currents at such depth.” 


Ripple marks and scour marks have been found on every seamount photo- 
graphed to date (Menard 1952, Heezen, ‘Tharp & Ewing 1959). Ripple marks have 
been observed on the crest of the Mid-Atlantic Ridge down to 2000 fathoms 
(Elmendorf & Heezen 1957). ‘They have also been observed on the continental 
slopes to similar depths. The deepest photographs of the North Atlantic abyssal 
floor had indicated probable scour marks in 3200 fathoms depth (Dietz 1952, 
Elmendorf & Heezen 1957). The deepest Atlantic photographs recently taken 
revealed rock bottom, at 4000 fathoms in the Romanche Trench portion of the 
Mid-Atlantic Rift Valley (Cousteau 1958). ‘These sharper topographic features 
are being scoured by ocean currents which remove the finer sediment and produce 
the ripple marks. The direction of these currents has not been determined, although 
a collection of oriented photographs might go far to solving the problem. It must 
be supposed that the sediment eroded from these peaks by the currents is deposited 
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in their lee. A study of these deposits may produce evidence of the direction of the 


deep currents and the nature of their variations in late Tertiary, Pleistocene and 
Recent time. 


Deep Currents. Wiist (1935, 1957) concluded that relatively strong bottom 
currents swept along the continental slope off eastern South America. His calcula- 
tions indicated a 5~20cm/s western boundary current in the South Atlantic. The 
same range of velocities was recently measured in the North Atlantic western 
boundary current through the use of neutral floats (Swallow & Worthington 1957). 

Ocean-bottom photographs of rock outcrops and ripple marks taken on the 
continental slope shows the eroding capability of this current (Northrop & Heezen 
1951, Heezen & others 1959). Dietrich (1957) and Cooper (1955) have described 
bolses of cold Norwegian sea water which flow through the Denmark Strait and 
descend to the depths of the Atlantic at velocities of 10 to 30cm/s. Strong currents 
have been calculated from the fine systematic observations made by J. B. Tait 
(1956) in the Shetland Channel. 

It has been suggested that the bolses of cold water flowing across the northern 
margins of the Atlantic may be related in some way to the turbidity currents which 
built or eroded the northwest Atlantic Mid-Ocean Canyon (Ewing & others 1953). 

Except where shallow-water benthonic organisms are found incorporated in the 
deposit, it may be difficult to distinguish in every case the effects of bottom scour 
from the deposits of turbidity currents. 

More must be learned of the causes and further effect of these strong bottom 
currents which have denuded the crest of a large part of the Mid-Atlantic Ridge 
to depths of 2000 fathoms (Heezen & others 1959). These currents wear so 
gradually that transatlantic submarine cables crossing the crest of the ridge often 
last 25 years or more before their armour wires are finally worn through. ‘Thus, 


this incessant motion must never become too violent or the cables would snap 
more frequently. 


Internal Tides. Internal tides have been observed during all oceanographic 
anchor stations made in deep water. Defant (1950) concluded that: 


“Internal tides are as common as tides at the surface of the ocean...’ ““The 


amplitude . . . and the water transport connected with them are much larger than 
those of normal surface tides.”’ 


In view of this it seems likely that internal tides could be responsible for the 
abyssal ripple marks observed in deep-sea photographs. 

Swallow (1955, 1957) reported eleven direct measurements of deep currents 
made in the eastern Atlantic with neutrally buoyant floats. Residual non-tidal 
components of 22cm/s down to 1cm/s were measured in depths ranging from 
340 to 2goom. A tidal current varying from 7-3 to 10°0cm/s was found super- 
imposed on a net current of 2-4cm/s in one observation at 800m depth. How- 
ever, in measurements made at greater depths the tidal component, if present, 
could not have exceeded 5 cm/s. Swallow suggests therefore: 


“Possibly the contribution from internal waves becomes less important at 
greater depths, so that the remaining tidal current may be near the theoretical 


value of about 3.cm/s derived on the assumption that the tidal current is uniform 
over the whole depth.” 
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If Swallow’s tentative conclusion is further confirmed it would appear that 
the deep ripples and bottom scour must be the work of currents related to the 


general circulation of the ocean rather than to the work of currents induced by the 
tides. 


Homogenized History. Evidence of deep-sea erosion, transportation and deposi- 
tion by bottom currents comes from the fine studies of the minute, fragile tests of 
silicious Radiolaria and diatoms found in mid-Pacific abyssal red clays. The 
long, red clay cores of the Albatross Expedition were chosen for detailed study 
since, based on estimates of sedimentation rates, they were thought to contain the 
undisturbed record of most of the late Tertiary. At the time Riedel (1957) began 
his investigation of the Radiolaria, he 


“expected that the assemblages in the tops of the cores could be used to 
determine the areal distribution of various species at the present day, while 
samples from lower levels of the cores would show, in an orderly manner, the 
development of present-day species through at least Pleistocene and Recent 
time.” 


However, the results of his study forced him to conclude: 


“that the situation is actually more complicated than was expected. Late 
Tertiary Radiolaria were found to be reworked into Recent sediments, and evi- 
dence of disconformities was found in the sediment sequences. Re-examination 
of some of the samples from the Challenger Expedition showed that ‘Tertiary 
sediments outcrop on the deep-sea floor. Erosion, and reworking of older Radio- 
laria into younger sediments, obscure phylogenetic developmental sequences in 


the long sediment cores . . .” 


The ‘‘epic poem” contained in these layers of sediment had, it seemed, been 
frequently reduced to a font of jumbled type as cryptic to read as a typesetter’s pie. 
Kolbe (1957a, 1957b) who found pre-Quaternary diatoms in several Pacific and 
Indian Ocean cores, concluded that their presence is “unimportant” and “‘acci- 
dental, although no plausible explanation of the occurrence can be given”. Riedel, 
however, “‘speculated that a more or less continuous agency (? current) was 
transporting Radiolaria on or near the sea floor” and that “the transporting 
agency has acted either continuously, or at such short intervals as to obscure 
intermittance”’. 

In the light of Riedel’s conclusions, it is not surprising that difficulty was 
encountered in the stratigraphic interpretation of the chemical analyses of these 
cores (Berrit & Rotschi 1956), and that, in general, the ionium and radium methods 
of dating of deep-sea sediments have been more often a failure than a success 
(St. Kroll 1955, Volchok & Kulp 1957). Probably the most important result of the 
later studies has been the conclusion that conditions of uniform, undisturbed, 
unreworked sedimentation required by these dating methods only rarely prevail 
on the deep-sea floor. 


Radiocarbon evidence. The top samples of deep-sea cores generally show an 
apparent radiocarbon age of the order of a few thousand years. Even when the 
“zero age” for the top of the core is extrapolated from several dates determined 
deeper in the core, ages are found averaging 2500-3000 years. It has been 
claimed that this discrepancy results from the action of bottom burrowers which 
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' churn the sediment bringing older material to the surface. Others believe that 
sedimentation underwent a recent drastic diminution. Isotopic fractionation 
during sedimentation, decomposition and diagenesis has also been suggested. 

The age determined for the fine fraction of a sample is generally greater than 
the age of the coarse fraction, and the age of the organic part is generally less than 
that of the carbonate portion. However, these generalizations frequently do not 
hold. 

Broecker, Ewing & Heezen (in press) attribute the finite zero age of top samples 
at least in part to the continuous erosion and redeposition of sediment. Sediment 
suspended in the water or wafted along the bottom is believed to be continuously 
reworked into the modern sediment. In addition, they believe that loss of sedi- 
ment during coring and churning by bottom burrowers are also significant 
contributory factors in many cases. However, the fact that similar finite zero ages 
are found in varve-like sediment rapidly deposited in anaerobic basins (Heezen & 
others 1958) indicates that these factors are insufficient as a general explanation. 


Sub-bottom Echoes. "Vhe character of the echo returned to the echo sounder 
from the deep-sea floor can be examined in great detail by use of the high-resolu- 
tion Precision Depth Recorder. 

The echo may be sharp and discrete, diffuse and prolonged, single, double or 
multiple. Over distances of a few hundred yards, the bottom may appear abso- 
lutely smooth or completely disrupted on a minute scale. 

Sub-bottom echoes are frequently observed in deep-sea echograms. In some 
cases the sediment horizon responsible for the sub-bottom echo has been identified 
by coring or drilling. Worzel (1959) has correlated a prominent sub-bottom echo in 
the east tropical Pacific with a 2—4 inch thick bed of volcanic ash. Heezen & others 
(1958) have correlated a prominent sub-bottom echo with the base of the Recent 
anaerobic layer in the Cariaco Trench. Gorsline & Emery (1959) have mapped 
turbidity-current beds in the San Pedro ‘Trough. In areas of particle by particle 
deposition where one or two prominent, 5-20-fathom deep sub-bottom echoes can 
be traced over a wide area they typically show deepening in the small depressions, 
thinning over gentle rises, and are totally absent on steeper inclines, indicating 
that due to gentle scour sedimentation has been greater in the depressions. 

The gentle but continuous wafting suggested by Riedel in explanation of the 
distribution of Radiolaria would seem to be a satisfactory explanation of these 
configurations. In other areas the simple sub-bottom echo, while still detectable, 
became intensely disturbed, possibly suggesting that the bed is intermittent or has 
been folded. 

The occurrence of a series of closely-spaced reflectors on an otherwise smooth 
bottom (Heezen & others 1959, Plate 5) may indicate deep-sea dunes (2 400 
fathoms) produced by strong bottom scour. Echogram profiles taken in depths of 
300 fathoms beneath the Gulf Stream were interpreted in a similar manner. 

Photographs of manganese nodules (Dietz 1952, Elmendorf & Heezen 1957) 
taken in great depths of the Atlantic Ocean floor, and more recently on the Easter 
Island Ridge of the southeast Pacific (Revelle 1958) tell us that sedimentation in 
these areas is virtually at a standstill. Pettersson has estimated their rate of 
concretionary growth at 1mm in 1000 years on the basis of measurements of 
radium. Regardless of their actual rate of growth, which may be one-tenth or 
one-hundredth of that amount (Goldberg 1958) it must be exceedingly slow for 
Tertiary sharks’ teeth often form their nuclei. Features in the sediment around 
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the nodules resemble scour marks. Thus the gentle scour and reworking of the 


abyssal red clay is supported both by the study of cores and by examination of the 
bottom photographs. 


Archaic Fauna. ‘That the bottom currents were not always so strong, so fast, 
or so cold, is shown by the frequent discovery of early Tertiary calcareous oozes 
beneath the red clay of the eastern Pacific (Arrhenius 1952). In several areas 
where now less than a few per cent of the carbonate falling onto the modern sea 
floor escapes solution, in the Eocene and Oligocene seas over seventy per cent of 
the carbonate remained. Oxygen-isotope measurements have been interpreted by 
Emiliani & Edwards (1953) to indicate cooling of 12 °C in the Pacific bottom water 
since Eocene time. However, regardless of the various assumptions of the isotope 
method, the existence of the high carbonate deposits itself is indicative of a com- 
parable change of temperature or rate of circulation. 

The conclusion of Menzies & Imbrie (1958) that the deep-sea fauna contains 
fewer rather than, as formerly held, a greater number of archaic types is probably 
the result of the apparently relatively wide variation in the rates and modes of 
circulation of the abyssal waters as contrasted to the almost uniform water condi- 
tions in the near surface environments. ‘The close uniformity of variables in the 
deep-sea with depth and geographical location had misled earlier investigators to 
assume that the deep-sea was an environment unvarying over long periods of time. 
But it now seems that the properties of the deep-sea most important to living 


organisms, namely, oxygen, temperature and nutrients, have varied widely 
throughout geologic time. 


Some problems 


Solution: The Big Lollipop. \n addition to mechanical processes of erosion, 
deposition and mixing by slumps, turbidity currents, deep-sea currents and tidal 
oscillations, the sediments of the deep-sea floor undergo chemical alteration and 
solution before, during and after burial. It is frequently incorrectly assumed that 
once buried the components of deep-sea sediments remain locked in the ocean 
depths, permanently removed from the bio-geochemical cycle. ‘The interstitial 
waters gradually forced from the sediment during compaction probably deliver 
some components back to the sea water (Emery & Rittenberg 1952). Volcanic 
glasses are devitrified and altered to clays and carbonates. Calcareous sediments 
are dissolved not only while sinking through the water, but after burial, since the 
mechanical processes active on the sea floor favour solution and alteration. 

In addition to currents and waves, the benthic fauna continually stirs the sedi- 
ments. Bottom photographs show abundant trails and tracks criss-crossing the 
sea floor. Often the bottom is studded with small cones a few inches across built 
by burrowing organisms (Heezen & others 1959). Sediment cores characteristically 
show burrow mottling. 

With abundant evidence of erosion, reworking and redeposition by bottom 
currents, we now must conclude that carbonate, silica and other components may 
re-enter the bio-geochemical cycle after millions of years residence beneath the 
deep-sea floor. 

Geochemical balance calculations generally encounter difficulty in regard to 
the high proportion of calcium in sedimentary rocks (Kuenen 1950). In explana- 
tion it has been suggested that no carbonate deposits could have formed in the 
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Palaeozoic ocean due to the lack of lime-secreting plankton; but Bramlette (1958) 
points out that cocoliths which make up the bulk of modern pelagic carbonate 
sediments, probably were present in Palaeozoic seas. However, the evidence of 
reworking, erosion and resolution of carbonate sediments long entombed beneath 
the sea floor may ease the problem since it need not be assumed that the calcium 
in deep-sea sediments is permanently lost from the geochemical cycle. 


**Moholes.”” One of the major problems of deep-sea sedimentation, indeed 
probably the most perplexing question in regard to the origin of the ocean floors 
is the thinness of the deep-sea sediments as determined by the seismic-refraction 
method. Over most of the Pacific (Raitt 1956, Hill 1957) the unconsolidated sedi- 
ment and sedimentary rock has a thickness of less than a half kilometre; in the 
Atlantic the average thickness of deep-sea sediments (Ewing & Ewing 1959) is 
only slightly greater. If sedimentation at rates at all comparable to those measured 
for post-glacial time by radioactive decay are applied to this thickness, we come to 
the startling but dubious conclusion that the oceans were formed in late Palaeozoic, 
late Mesozoic, or early Tertiary time; the exact date depending on our particular 
choice of average rates. This calculation allows for only pelagic sedimentation and 
takes no account of turbidity-current additions to the thickness. If we cannot 
assume that throughout most of geologic time the marginal trenches trapped the 
turbidity-current deposited sediments, the discrepancy becomes much greater. 
As an explanation of this apparent discrepancy, the solution of deposited sediment, 
the submergence of most of the present lands with attendant decreased erosion, 
and the lack of Palaeozoic plankton all have been suggested ; but the problem of the 
extreme thinness of the sediments has become nonetheless evasive and intriguing. 
The answers to these problems may be provided by core samples obtained during 
the drilling in the deep-sea floor of the exploratory holes to the Mohorovici¢ 
discontinuity (Gaskell 1958). However, in view of the many phenomena which can 
continue to jumble the record, great care must be exercised in the choice of 
drilling sites to better our chances of obtaining a readable record. 


Conclusion 


We now can view the vast abyss with somewhat less mystery and more real 
understanding. The picture of an unremitting snowfall of sediment which softly 
casts a thin veil over the rugged, grand and imposing relief of the deep-sea must 
be tempered with occasional, sometimes frequent, turbidity flows, which, like 
snow avalanches in the Alps, sweep down and smother the relief of the lower 
slopes. Just as winds striking the mountains drive the snow from the peaks, so 
the ocean currents sweep the sediment from the higher and bolder submarine 
mountains. 

Just as wind and water shape the land so do the submarine currents of water 
and mud, by erosion, by solution and by scour shape the ocean floor. To be 
sure, erosion dominates on the land, as does deposition in the deep sea. 
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Reports of Discussions 


Geophysical Discussion of the Royal Astronomical Society 


A discussion on The Characteristics of the Antarctic Ice Sheet was held in 
the rooms of the Society on 1959 February 27. The meeting was arranged jointly 
by the Royal Astronomical Society and the British Glaciological Society under 
the chairmanship of Dr J. F. Nye. 

After welcoming and introducing the other speakers Dr Nye gave an outline 
of the theoretical aspects of the problems presented by large ice sheets. The 
mechanical behaviour of an ice sheet is governed by two opposing forces: the 
weight, which tends to squeeze the ice out sideways, and the resistance to outward 
flow provided by the rock bed. These two forces are approximately balanced. 
Two substantial sources of heat affect the temperature distribution within the 
ice. Heat flows into the bottom of the ice sheet from the earth itself and heat is 
generated at the bottom of the ice sheet by friction between the ice and the rock 
bed over which it moves. If the minor effects of transverse convection are ignored 
both sources of heat tend to raise the temperature of the lowest layers. As the 
plasticity of ice is sensitive to small changes in temperature, it follows that the 
plasticity of the ice sheet is greatest near the bottom. The shearing stress in the 
ice increases linearly with the depth and hence the shear distortion takes place 
predominantly in the lowest layer and analysis has shown that the movement may 
equally well be regarded as slipping on the bed or as very intense shear in the 
lowest layer. The two principal features of the motion are thus a very slow 
squeezing down partly or entirely compensated by the accumulation of new snow 
and giving a strain rate of about 10~4/yr, and a faster outward movement giving 
a shear rate greater than 10~2/yr in the bottom layer. In a steady state the total 
outward flow of ice is just sufficient to balance the accumulation of new snow. A 
simple model shows that the steady state thickness of the ice sheet is proportional 
to only the 5th or 6th root of the accumulation and so the height of the ice sheet 
will be very insensitive to even long term fluctuations in the rate of accumulation, 
and calculations of the form of an ice sheet can be made with fair accuracy without 
a knowledge of the accumulation. Field measurements on a large ice sheet should 
aim at elucidating the topography, ice thickness, and rate of accumulation. 

The second speaker was Ingénieur Hydrographe B. Imbert who is in overall 
charge of the French I.G.Y. Expeditions to the Antarctic and who himself led the 
second French I.G.Y. Expedition to Terre Adélie. M. Imbert presented the 
results obtained along a traverse running 300 miles inland from Base Dumont 
d’ Urville to a southernmost point south of Base Charcot at an altitude of 2 509 m. 
During the southern summer 1957-58 seismic reflection results were obtained at 
twenty-five stations along this route and during the next summer seventy gravity 
stations were occupied and glaciological work was carried out. 

A twelve channel recording camera and filter with pass bands of 5—100c/s 
were employed in the seismic work. From the coast to an altitude of 480m 
reflections were obtained with shots fired on the surface but farther inland it was 
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necessary to bury the charges and at the southernmost point, where the ice was 
2990 m thick, a charge of 1kg was used at a depth of 18m. Larger charges gave 
no improvement. Difficulties were experienced from the cold surface layers of 
the snow which tend to spoil transmission and from persistently strong winds 
averaging 30 knots. It was impossible to record in winds greater than 38 knots. 
In spite of this the typical reflection records passed round by M. Imbert were of 
a very high quality. The results of the seismic work showed that the undulating 
bedrock is largely below sea-level along the line of the traverse. The average 
depth is 40m below sea-level but a maximum depression of 600 m below sea-level 
was recorded. M. Imbert suggested that a broad valley observed near Charcot is 
occupied by a glacier flowing to the north-east because instruments levelled at 
Charcot gradually became tilted in this direction. 

In establishing the topography of the ice sheet for the gravity and seismic 
traverses heights were estimated barometrically but the atmospheric pressure was 
checked independently by radio-sonde methods. The estimated accuracy of the 
heights obtained is +15-20m. ‘The results of the gravity traverse are not yet 
available. The profile of the ice sheet along the line of the traverse is given by the 


ellipse 
h (550—x\2 
2°68 550 
where x is the horizontal distance in kilometres inland from Dumont d’Urville 
and h is the height in metres of the ice surface above sea-level. The r.m.s. error in 
his 25m. 

Glaciological work was carried out by the French expedition. Fifty-seven 
stakes were set up on the traverse of 483 km and were observed during the two 
subsequent summers. More detailed observations of accumulation were made at 
Charcot. The results showed that the rate of accumulation falls off from a value of 
36cm of water equivalent per year near the coast to 10cm per year at the southern- 
most point. The observed surface densities lay between 0°38 and 0-48 g/cm’. The 
mean temperature of the ice at Charcot was — 38°C at a depth of 16m. 

In answer to Dr T. F. Gaskell the speaker stated that multiple explosion shots 
fired above the snow surface did not improve the records sufficiently to justify 
the increased expenditure of explosive. The Chairman enquired what was the 
magnitude of the tilting observed at Charcot and M. Imbert estimated that it 
was a few minutes of arc per fortnight. Dr H. Lister noted that the rate of accumu- 
lation observed at Charcot was high and enquired whether the station was in a 
topographic trough. M. Imbert said that the station was not in a trough but 
agreed that accumulation rates were greatest in troughs. 

The third speaker, Mr J. G. D. Pratt, described the seismic and gravity 
traverse carried out by the Commonwealth Trans-Antarctic Expedition. Mr Pratt 
said that, like M. Imbert, he had found a range of 5—100 c/s adequate for the seismic 
recording equipment but that it had been necessary to fit a heater into the camera. 
Lack of time had prevented him from burying the shots more than 12m below 
the surface and as this produced no great improvement the shots were sometimes 
fired on the surface. A slide showed records obtained at the South Pole where 
there is a discrepancy between Mr Pratt’s estimate of the depth to the bedrock 
and the estimate based upon a single American observation. Mr Pratt pointed out 
the reflection which he ascribed to the bedrock and the subsequent reflection 
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which the American observers had chosen. He said that his own estimate of the 
depth, which agrees with the results of neighbouring seismic stations and with the 
gravity profile, should be accepted in preference to the very much greater American 
estimate. Mr Pratt showed a profile of the ice sheet from South Ice on the Weddell 
Sea side of the Pole to the Plateau Depot on the Ross Sea side. The results 
obtained by the Trans-Antarctic Expedition showed that along the line of their 
route the rock bed was almost everywhere above sea-level. On the Weddell Sea 
side there are two depressions comparable with Lake Baikal, for example, but on 
the Ross Sea side of the Pole the bedrock is substantially higher, due probably to 
the fact that the route of the Expedition was parallel with the known mountain 
range fringing the Ross Sea and probably lay over the buried foothills of the 
range. 

The results of a refraction line shot at South Ice were used to find the varia- 
tion in the velocity of compressional waves with depth in the ice. It was found 
that the velocity increased in the surface layers from 0-5 to 4m/s at a depth of 
80-g0m _ Below this the velocity was uniform. The observed rate of increase in 
velocity in the surface layers was used to calculate an accurate value for the cor- 
rection due to the depth of shot and it was shown that this correction could have 
been made with sufficient accuracy by the usual assumption of a linear increase 
in velocity with depth in the surface layer. Although the snow surface was undu- 
lating no correlation was found between the observed arrival times of the refracted 
waves and the topography. This was interpreted as meaning that the valleys are 
underlain by softer snow in which the lower velocity counterbalances the effect of 
the shorter path which the sound has to follow. 

Shackleton Base was built on an ice shelf. Reflection shots were fired near the 
Base to try to determine the thickness of the floating ice. Reflections were obtained 
from the sea bed but no reflections were obtained from the bottom of the ice 
shelf. However the travel times for a ray refracted at the ice—-water interface and 
reflected at the sea bed can be expressed as a power series in terms of the hori- 
zontal distance from shot to detector and Mr Pratt has shown that a unique solution 
of the thickness of the ice shelf can be obtained in this way. 

The results of the gravity measurements were used to confirm the general 
correctness of the seismic profile. Altitudes were measured barometrically and 
errors from this source amount to +25-somgal. The estimated errors due to 
drift are about 5mgal. Bouguer anomalies have been calculated taking the ice 
thicknesses indicated by the seismic profile and assuming a rock density of 
2-1g/cm3. The fact that there are no outstanding anomalies is taken to indicate 
that the seismic profile is substantially correct. However there is a mean anomaly 
of — 200mgal which cannot be explained away by the uncertainties of the measure- 
ments and suggests that Antarctica, like Scandinavia, has not yet recovered from 
the depression caused by the burden of a much larger ice sheet which covered 
the continent in the past. 

The last of the four speakers was Dr H. Lister, who described the glaciological 
results obtained by the Trans-Antarctic Expedition. In order to determine the 
mass balance of the ice sheet as a whole it is necessary to know how the average 
annual rate of accumulation varies from the centre of the continent to the peri- 
phery. The average annual rate of accumulation at the coast was estimated as 
equivalent to about 30cm of water per year by observing the layering in the 
shelf ice cliff. In order to observe the accumulation during the year at Shackleton 
Base and at ‘‘Southice”’ strings and graduated pegs were fixed to the snow surface 
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between stakes. This was done every fortnight. At the end of the winter the strings 
were dug out and the fortnightly accumulation measured. In order to find the 
rate of accumulation in past years, observations were made of the layering in pits 
and boreholes and the hardness of the snow layer was measured by means of the 
Rammsonde device. The annual layering is marked by variations in grain size 
hardness, and density, but it is not easy to be certain of the identity of the individual, 
layers. In pits and boreholes at ‘“‘Southice” Dr Lister was able to recognize the 
annual layering back to the year 1895 with some confidence and the mean rate of 
accumulation obtained was about 1ocm/yr. On the trans-continental journey 10 
boreholes penetrated through twenty annual layers and 40 more penetrated five 
layers. Zones of hard sastrugi and areas of broad undulations were encountered on 
the Polar Plateau. The mean rate of accumulation found on the plateau was 8-10 
cm/yr but as much of the route of the expedition lay in a topographic trough this 
figure is probably a little greater than average. In the coastal mountains and on the 
shelf ice the observed rate was 16cm/yr rising to 26cm/yr near the sea. ‘The low 
accumulation observed as Shackleton—10cm/yr—was explained as being due to 
loss of drift snow over the nearby ice cliff. Making due allowance for the shape of 
the continent Dr Lister calculated a mean annual rate of accumulation of 13-1 
g/cm? yr for the whole Antarctic ice sheet. 

Photographs and observations were made to determine the distribution of drift 
snow above the surface and led to a figure of 0-1 g/cm? yr for the loss from the 
whole continent due to drift blown out to sea. Assuming a mean flow rate of 
20o0m/yr for the edge of the ice sheet Dr Lister arrived at a mass balance for the 
continent of 13-1g/cm? yr accumulation and 10-1 g/cm?yr total loss due to evapora- 
tion (1-0g), deflation (0-1g), calving of icebergs (4-og), and melting, largely at the 
ice/sea interface (5-0g). He said that there was probably a closer approximation 
to balance than these figures would suggest. 

Mr Pratt pointed out that it was curious that the tracks of drift particles shown 
on the photographs were all parallel, indicating that all the particles were at the 
same point on their trajectories. Dr H. H. Lamb of the Meteorological Office 
asked what were the estimated errors of the points on the curve showing annual 
snow accumulation at “Southice” since 1898. Dr Lister said that the standard 
deviation of a single observation was 34g but that the actual error probably amoun- 
ted to 50 per cent. 

In the discussion which followed the papers Dr T. F. Gaskell asked how it 
might be possible to reconcile Dr Lister’s results which suggest net accumula- 
tion with Mr Pratt’s gravity results which suggest net loss. Mr Pratt said that the 
gravitational results indicate a long term retreat, whereas the glaciological results 
apply to a short period of only a few years and may indicate a minor fluctuation. 
The Chairman suggested that the rate of flow of the ice sheet should be taken as 
greater than 200m/yr and Professor J. H. Zumberge said that he had measured a 
flow rate of about 10 metres per day on the Ross Ice Shelf. Dr Lister said that it 
would not be sufficient to double the estimate of 200 m/yr which he had taken but 
that the high flow rate of about 1000m/yr suggested by Professor Zumberge 
would be enough to remove the discrepancy between accumulation and loss which 
his figures had indicated. Professor Zumberge asked why the low figure of 
1 g/cm? yr had been taken for evaporation. Dr Lister replied that he had actually 
observed rime accumulation rather than evaporation. He thought that there is 
insufficient heat available for much sublimation and that 1 g/cm yr was a generous 
estimate of the loss from the continent as a whole due to evaporation. Professor 
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Zumberge mentioned measurements of summer meltwater on land, made in the 
McMurdo Sound area, but Dr Lister thought that this melting merely balanced 
the excessive snowfall which occurred near the sea. 

The Chairman asked what was the significance and history of the large-scale 
surface undulations observed on the polar plateau. Dr Lister said that in the walls 
of the pit dug at ‘‘Southice”’ a regular fluctuation in grain size had been observed 
and that this had been ascribed to the migration of large sastrugi with a frequency 
of about one per 14 metres of accumulation. This was comparable with the 
surface undulations which have a wavelength of about 12km and a mean amplitude 
of 20 to 25m. The Chairman pointed out that these undulations may have a pro- 
found effect on the motion of the ice sheet. They represent fluctuations in slope 
comparable with the average slope itself and extending over distances comparable 
with the thickness of the ice sheet. Since the shear stress on the bed is propor- 
tional to the surface slope and on Weertman’s theory the sliding velocity is propor- 
tional to the square of the stress, the undulations may result in a system of random 
velocities, perhaps greater than the mean velocity, superimposed on the smoothly 
varying outward velocity. This might invalidate measurements of the velocity of 
flow made in the field. Mr G. R. Elliston asked whether it would be possible to 
measure the grain size in a pit at the crest of an undulation in order to see how 
long the crest had persisted. Dr Lister replied that the random variation would 
make this impossible. Replying to Dr Zumberge he said that the undulations were 
asymmetric in profile on the Ross Sea side of the Plateau. Mr Hollingworth asked 
why the undulations did not smooth themselves out and what was the mechanism 
of their motion. Replying, the speakers agreed that bad sastrugi were found on 
the flanks of the undulations and that these might interfere with the flow of air and 
be responsible for the soft snow and greater accumulation found in the troughs. 

The Chairman thanked the speakers and the contributors to the discussion 
and adjourned the meeting at 18h 20min. 


The Geological and Geophysical Results of the cruises of 
R.R.S. Discovery II during 1958 May—July 


A symposium concerned with the geophysical and geological investigations of the 
ocean floor during 1958 from R.R.S. Discovery J] was held at the Department of 
Geodesy and Geophysics at Cambridge on March 4th. The expedition, which 
was a joint venture of the Department of Geodesy and Geophysics and the National 
Institute of Oceanography, lasted approximately three months and during this 
time investigations of the ocean floor were made in three separate areas. The first 
was in the neighbourhood of a small sea-knoll in latitude 41°15’ N and longitude 
14°30’ W, near the western boundary of the Iberian Abyssal Plain which lies at a 
depth of about 5 300m. The main object of the cruise in this neighbourhood was 
the measurement of deep ocean currents; this was not a full-time occupation and 
allowed time and opportunity for detailed investigations of the topography and 
magnetic field of the abyssal plain and neighbouring sea mounts in an area approxi- 
mately 60 miles from north to south by 100 miles from east to west. Seismic 
refraction shooting experiments, heat flow measurements, coring, dredging and 
bottom photography were also successfully undertaken. 

The dredging over the sea-knoll and over the sea-mounts to the westward 
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provided a large collection of rocks which were not erratics and which, there was 
good reason to suppose, were in situ. Bottom photography had shown the existence 
of outcropping rock and the samples from any one haul were largely of the same 
type. The rocks from the sea-knoll were rotted basaltic lavas which were more or 
less manganese coated; this coating in some cases was entering cracks in the rock 
and causing the rocks to fracture. Weathering of this type has not hitherto been 
observed. Only small samples of unaltered glass were obtained but these will be 
adequate for chemical analysis; sections of these rocks have been made and are at 
present being studied. 

The lavas from the sea-mounts were less altered and of a different type from 
the others. Petrographical work on these rocks is now being undertaken. Measure- 
ments of the radioactivity, the magnetic susceptibility and the elastic properties 
are also to be made in the future. 

The seismic results in the area suggested that the Mohoroviéié discontinuity 
was at a relatively small depth. The combined results from three stations showed 
that on the abyssal plain a thin layer (about 200m) of sediments covered a iayer 
which, with an assumed velocity of 4km/s, was 2-7km in thickness. This overlay 
the typical oceanic “basaltic” layer of velocity 6-8km/s and thickness 1-3 km. 
The Mohorovitié discontinuity lay at a depth of 9-5 km below sea-level and the 
velocity below was 8-1 km/s. 

Various minor features were observed to stick up through the abyssal plain 
and the sediments collecting on their slopes were fine clay and foraminiferal tests. 
This material was markedly different from that in the cores from the abyssal plain 
which were loaded with coarse quartz sand and fragments of metamorphic rocks. 

The total magnetic force contours of the area showed a restricted correlation 
with the topography. The central sea-knoll had a clear anomaly associated with 
its northern end, and the sea-mounts to the west likewise caused greatly increased 
magnetic roughness. Other features on the abyssal plain, however, had no anoma- 
lies associated with them and conversely there were a number of magnetic anoma- 
lies dissociated from topographic features. From these measurements it was not 
possible to decide whether the anomalies were caused by rocks immediately below 
the sediments or by deeper rocks. 

The precision echo sounder developed by the National Institute of Oceano- 
graphy, and which was first used on this expedition, allowed great detail of the 
fine structure of the abyssal plain to be studied. It was found to slope upwards 
towards its boundaries to the north and east with slopes of a few parts in 10000. 
This had been deduced from earlier work. 

The second area of investigation lay in the region between the Biscay and the 
Iberian abyssal plains. This area was selected since it appeared from the slopes 
of the plains that there must be a channel or channels interconnecting the two 
plains through the saddle separating them and that sediment was transported from 
the Biscay plain to the Iberian plain. The depth difference between these plains 
is about 200m. Two channels were located on the Biscay side of the saddle; 
these joined into a single channel approximately halfway between the plains and 
then bifurcated into two at the Iberian plain end. The total distance between the 
plains is about 30 miles. The gradients in the beds of the channels are high at their 
top ends, slopes of 1 in 10 existing near the spill-over from the Biscay plain. 
Lower down, the channels have relatively small gradients and near the region 
where they join the gradient is in fact reversed. It is postulated that turbidity 
currents from the Biscay plain acquire sufficient kinetic energy when passing down 
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the steep slope to sweep over this region of slightly reversed gradient before travel- 
ling far out into the Iberian plain. Core samples from the bed of the valley were 
of coarse material in comparison with that obtained from the areas elevated above 
the abyssal plains. 

The third area of investigation was centred at a position 42}°N, 114°W 
and has provisionally been named Galicia Bank. It is an elevated plateau of depth 
about 7oom separated from Cape Finisterre by a channel of depth greater than 
z200om. ‘To the north lies the Biscay abyssal plain and to the south-west the 
Iberian plain. The plateau which is about 30 miles in length from north to south 
and 20 miles in width from east to west is bounded by steep escarpments to the 
east and to the west. Magnetic measurements show that, compared with measure- 
ments made over sea-mounts further from the continental boundary, the rocks 
must have a relatively very low magnetic susceptibility. In this respect the rocks 
appear to resemble rocks of the continental shelf more closely than those of the 
deep ocean. Seismic measurements indicate that the plateau is covered with a thin 
layer of unconsolidated material overlying rocks with velocity 4-8km/s. This 
layer is approximately 4 km in thickness and lies above a layer of velocity 7-1 km/s. 
From these velocities it is not possible to identify the rocks; they are not, however, 
inconsistent with what might be found in a column of continental rocks. Dredging 
in the area produced a collection of stones which, to date, have not been closely 
inspected. ‘They included a considerable number of metamorphic rocks, some 
basic lavas and some limestones. It is not yet certain that these rocks are in situ; 
it is possible that they are erratics. 

During the expedition numerous unoriented cores up to 7m in length were 
obtained for palaeomagnetic measurements. ‘The measurements are not yet com- 
plete but it appears from the preliminary results obtained from three cores taken 
from the flanks of a small sea-mount in the Iberian abyssal plain that they are 
magnetically stable. Correlation between the magnetic properties from core to 
core is not immediately obvious; future measurements may nonetheless show that 
it exists. ‘The mean value of the inclination of the magnetic vectors in these cores 
is not appreciably different from what would be expected if the present inclination 
and the associated variation of direction of the Earth’s magnetic field had persisted 
throughout the time of deposition of the core. The intensity of permanent mag- 
netization of the core material is directly associated with the water content which 
is itself associated with the quantity of clay minerals present. It appears that the 
increase of the proportion of clay results in an increase in water content. Extrac- 
tion of the magnetic minerals and a closer study of the lithology of the cores will 
take place in the future. 

During the expedition two refuelling visits were made in Lisbon and a very 
large magnetic anomaly was located in shallow water a few miles outside the mouth 
of the Tagus. A survey of the area showed that the anomaly could be produced 
by a dipole buried at a depth of approximately 2km below the sea floor. This 
dipole could be caused by a spherical mass of magnetic rocks of maximum radius 
2km centred at a depth of 2km. The minimum susceptibility of such a mass of 
rocks must be exceptionally high. The anomaly was associated with magnetic 
compass deviations of about + 4°. 

Those participating in the Symposium were: Dr M. N. Hill, Sir Edward 
Bullard, Mr T. D. Allan, Mr M. J. Keen and Mr D. H. Matthews, from the 


Department of Geodesy and Geophysics, and Dr A. S. Laughton from the National 
Institute of Oceanography. 
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KOJIEBAHHA 3EM.IH CBASAHHBIE C OTHM BOLPOCbI 
X. Tarxaymu 


yupyrux B KOTOpPLIX MeHAIOTCH MAOTHOCTH KOHCTAaHTLI, a 

Ckopoctn Bemau no Byaseny (cm. dur. 1). 
pacnpeyereHHA Ha dur. 2. 


MCKA#KEHHA, BHOCHMBIE CEACMOPPA®MAMH B 3EM.TETPAICEHHA 
MMITY.IBCHOPO BIA 


M. Jlanducman, H. Camo, u. M. lOune 


NHOCTOAHHLIMM B CHPHAIbE MMITYJIbCHOPO TeopeTH4eCKH Tak M 
HEX 060m ceiicmorpade 1A HenpephBHOTO 


veTBeptTore THHa Pynre-Kytra. Hamenenne CymlecTBeH- 

1 1938 Bp Hosoit Teunee ¢ ux we sanuchw Ha Aedopmorpade. 
PeayabtaTh Takoro H3y4eHHA C CaTo 0 TOM, 4TO rpaduK rpynnoBolt 
CKOpOCTH BO.IH COBEPIICHHO NPAMOJIMHEECH B LIMPOKOM NepHOROB B HECKONbKO 
coten cekyHa. Hekaskenne, BHOCHMOe npH6opaMN, KOTOpHX OLicTpo BO3pacTaeT 
FOCTATOYHO AIA TOPO, YTOOM MACKHpoBalb YO 
KoTOpan OLITL CBAZAHA C TAKOrO Poa 


ONPE]IEJIEHHE TEMMEPATYP B HEHAPYINIEHHbIX HO 
HAB-TIOJEHHAM B PIYBORHX PA3SBEJIOUHbIX 


I. P. Kynep u C. 


B npouecce OypeHuA CKBaxKMHEI TeMMepaTypa ee HapyulaeTcA 3a 
Gypoporo pacTBopa M OcTaHOBKH GypeHHA AKMAKOCTh B 
M kK Heil npHoOpeTaWwT TemnepaTypy 
TeMMepaTypHOro paBHOBeCHA CKOpee BCero nponcxoyuT y 38a60n 
CKBaKMHH, rye Hapywenne Ha OCHOBaHMM 9TOTO MOKHO 
CHMMeTPHH C YCAOBHAMH, TeNAOBLIMH NapaMeTpaMM 
it 

Paspa6oTana HOBaA M HeOporad annapatypa, TepMACTOpa 
H3MepeHHii, KoTOpan ycnemHo onpoSopana rayOun 1060 B cKBaxKuHe 
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Xontou Ilya patione Kannok Ueiic. Onncano M MaMepeHHbie TeM- 
hepaTyph M1acTOB, BCKpbITHX CKBaxKuHOM Xontou Ilya. cepna 
Temnepatypht Gyponoro pacTBopa Ha raiyOunax, B 
TeYeHHe Hel CHYCTA OKOHYAHMA M ITO 
OTH HE MOTYT Kak O TeMNepaType HeHapy- 
11acTOB. 


BOJIH B CPEJIE C MAJIBIMH TIOTEPHMIL 
Poce 


(IA MH BOTH B MaTepHadax, 
notepamu. pesyabTaTh H KOTOpHe NOsBOAAWT 
K JIKOOOMY, peadIbHOMY ABYX Tpex 
cKopoctel. 


KPATROE COOBLITEHHE 


B HHTEPPASIE CTOKCA 


4. Brammauapdoeu 


sonannn (popMy.tht Crokca BbIACHeHHA 
Semi HOBOLO BHA HOBepXHOCTHOM aHOMAIHM BHAA g—y,, 


O TTPOPPECCE B 


JMUHAMUYECKHE HPOWECChE PIYBOROBOJLHOTO OCAJTROOBPASOBAHHA 
SPOBHA, BUAKEHHE HAHOCOB H HA PTY BORKOBOJLHbIX 
YYACTKAX MOPH 


B. C. Xusen 


BOHOM Kak Ha CHOKOTHOe, 
B BpeMA 9TO nepecmotpeHo. mpoucxoqaume na 
CTBMeM MYTbEBLIX MOTOKOB, BHYTPCHHHX PIYOMHHEIX TeYeHH, pacrBopeHHA 
Ha OTHOCHTeIbHOe BHAYeEHHe TeX HIM HHBIX mponeccos. pasBMTHe HOBBIX 


M ONMCHIBAWTCH BAPHAKTH TOWeK BpeHHA Ha 
aTOT BOTIPOC. 
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